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ABSTRACT 

This  paper  presents  the  essential  features  of  Horowitz’s  transparency  of  QFT  and  his  QFT  technique  that 
exemplify  the  concept  of  “Bridging  the  Gap  ”  which  are  the  essential  aspects  of  the  QFT  control  system 
design  process. 
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PREFACE 

The  intent  of  this  paper  is  (1)  to  provide  a  basic  understanding  of  QFT  with  a  minimum  amount  of 
mathematics;  (2)  introduce  the  reader  to  the  QFT  design  technique  for  MISO  and  MIMO  analog 
and  discrete  control  systems;  (3)  to  provide  a  design  example,  and  to  illustrate  the  “Bridging  the  Gap.” 
Note  that  all  figures  in  this  paper  are  from  Houpis,  et  al.,  (1999). 


1.0  QFT  DESIGN  TECHNIQUE  FUNDAMENTALS 

QFT  is  a  frequency  domain  technique  for  the  design  of  robust  multivariable  control  systems  containing 
structured  and  unstructured  parametric  uncertainties.  This  paper  deals  with  the  former  case.  It  is  a 
powerful  multivariable  design  method  for  plants  with  structured  parametric  uncertainty  [see  Horowitz, 
(1991)  and  Horowitz,  et  al.,  (1973)],  for  systems  with  control  effector  failures  and  man-in-the-loop  flight 
control  system  designs.  QFT  also  accounts  up  front  the  plant  uncertainty  (LTI  plant  models),  performance 
specifications  (P.S.),  and  design  limitations.  It  also  has  the  feature  of  limiting  over  design  and  achieves 
reasonably  low  loop  gains,  i.e.,  avoids  or  minimizes  sensor  noise  amplification,  saturation,  and  high 
frequency  uncertainties.  QFT  technique  utilizes  a  simple  unity-feedback  control  system  structure  to 
achieve  a  desired  robust  system  design. 

1.1  Introduction  to  QFT 

The  material  in  this  section  presents  an  overview  of  QFT  in  order  to  prepare  the  reader  for  the  more 
detailed  presentation  of  QFT  in  this  paper. 

Why  Feedback?  —  It  is  necessary  that  the  feedback  control  system  design  technique  to  address  at  the 
onset  all  known  plant  variations,  to  incoiporate  the  information  on  the  desired  output  tolerances,  and  to 
maintain  a  reasonably  low  loop  gain  in  order  to  minimize  sensor  noise,  etc. 
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ORGANIZATION 


Consider  the  MISO  Control  System  of  Fig.  1(a)  —  In  this  figure  the  plant  uncertainty  set,  IF,  represents 
a  plant  with  variable  parameters,  the  desired  input  signal,  r(t),  that  is  to  be  tracked,  and  the  external 
disturbance  signal,  d{t),  to  be  attenuated  and  to  have  minimal  effect  on  the  output,  y{t). 

Tracking  and  Disturbance  Control  Ratios  (Open-Loop  System)  —  The  tracking  and  disturbance 
control  ratios  of  Fig.  1(a),  based  upon  the  nominal  plant  P0  e  IP,  are,  respectively: 

Y(s) 

Td(s)  = - =  P0(s)  (la) 

D(s) 

Y(s) 

Tr(s)  = -  =  P0(s)  (lb) 

R(s) 


Sensitivity  Functions  —  The  sensitivity  functions  of  Fig.  1(a)  for  the  two  cases: 
Yr(s)  I  dft)  =  o  &  Yd(s)  I  r(t)  =  o  are  identical;  i.e.: 


(s)  =  S 


ydW 

Vs) 


=  1 


(2) 


Tracking  and  Disturbance  Control  Ratios  (Closed-Loop  System)  —  The  tracking  and  disturbance 
control  ratios,  where  L0=  GP0,  the  nominal  loop  transmission  function,  of  Fig.  1(b)  are,  respectively: 


T  GP0  L0 
Tr  = - = - 

1  +  GP0  1  +  L0 


(3) 


Td 


Po  Po 

1  +  GPo  1  +  Lo 


(4) 


Sensitivity  Function  for  Compensated  System  —  The  sensitivity  function  of  Fig.  1(b)  for  these  two 
cases  are  also  identical;  i.e.: 


Yr(s) 

'PqOO 


(s)  =  S 


YDW 

PoW 


(s)  = 


1  1 
1  +  GP0  ~  1  +  L0 


(5) 


Comparing  Eq.  (5)  with  Eq.  (2)  illustrates:  (a)  the  effect  of  changes  of  the  uncertainty  set  T(.sj  upon  the 
output  of  the  closed-loop  control  system  is  reduced  by  the  factor  1/[1  +  GP0]  compared  to  the  open-loop 
control  system,  and  (b)  that  this  reduction  is  an  important  reason  why  feedback  systems  are  used. 


R 


Y 


(a) 
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(b) 

Figure  1:  Control  Systems:  (a)  Open-Loop  Compensated  System; 
(b)  Closed-Loop  Compensated  System. 
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Horowitz  —  Horowitz  has  stressed  that  the  best  robust  control  system  design  is  achieved  by  working  with 
L0  and  not  with  the  sensitivity  function  S. 

Reasons  for  the  choice  of  L0  —  The  reasons  for  the  choice  of  L0  are:  (1)  the  sensitivity  function  is  very 
sensitive  to  the  cost  of  feedback;  (2)  a  practical  optimum  design  requires  working  to  the  limits  of  the 
system’s  performance  specifications  (P.S.);  and  (3)  the  order  of  the  compensator  (controller)  G  can  be 
minimized  by  incoiporating  the  nominal  plant  P0  into  L0. 

Basic  Explanation  of  Structured  Parametric  Uncertainty  —  The  QFT  design  objective  is  to  design 
and  implement  a  robust  control  for  a  system  with  structured  parametric  uncertainty  that  satisfies  the 
desired  P.S. 

Parametric  Uncertainty  -  A  laboratory  experiment  involves  hooking  up  the  d-c  shunt  motor  as  shown  in 
Fig.  2  by  students.  When  they  entered  the  laboratory  on  a  cold  January  Monday  morning  to  perform  this 
experiment  the  weekend  room  temperature  was  50°F.  It  was  immediately  reset  to  70°F.  The  students  then 
hooked-up  the  motor  and  set  the  field  rheostat  to  yield  a  speed  co  =  1200  r.p.m.  Having  accomplished  this 
they  decided  to  take  a  one  hour  break  in  order  to  allow  the  room  to  reach  the  desired  temperature  of  70°. 
Upon  returning  they  found  that  the  speed  of  the  motor  was  now  1250  r.p.m.  with  no  adjustments  to  the 
applied  voltage  or  of  the  field  rheostat  Rf. 


Nonlinear  Plant 


Why  the  change  in  speed?  -  Due  to  the  heating  of  the  d-c  shunt  field  by  the  field  current  if  and  the 
environment,  the  value  of  Rf  increased.  This  increase  in-tum  decreased  the  value  of  if  and  in-tum  resulted 
in  an  increase  in  the  speed  since  it  is  inversely  proportional  to  if  assuming  Vf  is  constant.  Therefore,  during 
the  operation  of  the  motor,  the  parameter  Rf  can  vary  anywhere  within  the  range  R|rnMn  <  Rf  <  Rfmax  due  to 
the  variable  environmental  temperature  and  the  field  current.  As  a  consequence,  there  is  uncertainty  as  to 
what  the  actual  value  of  the  parameter  Rf  will  be  at  the  instant  a  command  is  given  to  the  system. 
Thus,  the  parametric  uncertainty  is  structured  because  the  range  of  the  variation  of  Rf  is  known  and  its 
effects  on  the  relationship  between  Vf  and  co  can  be  modeled. 


A  Simple  Mathematical  Description  of  a  d-c  Servo  Motor  —  The  transfer  function  of  a  d-c  servo  motor 
is: 


Pi(s)  = 


©m(s) 
Vf  (s) 


Ka 

s(s  +  a) 


(6) 


where  the  parameters  K  and  a  vary,  i.e.,  K  e  (Kmin,  Kmax)  and  a  e  (amin,  amax)  over  the  region  of  operation. 
For  this  example  Kmin  =  1  and  Kmax  =10  and  amin  =  1  and  amax  =  10.  The  plant  parameter  variations  are 
described  by  Fig.  3  where  the  shaded  region  represents  the  region  of  structured  parametric  uncertainty 
(region  of  plant  uncertainty). 


RTO-EN-SCI-142 


1  -3 


Quantitative  Feedback  Technique  (QFT):  Bridging  the  Gap 


ORGANIZATION 


Figure  3:  Region  of  Plant  Parameter  Uncertainty. 


The  motor  can  be  represented  by  J  =  6  LTI  transfer  functions  (v  =  1,2,... ,J)  at  the  points  indicated  on  the 
figure.  The  Bode  plots  for  these  6  LTI  plants  are  shown  in  Fig.  4  which  graphically  illustrate  the  structured 
parametric  uncertainty  in  magnitude  (dB)  and  in  phase,  i.e.,  for  a  given  00  =  ©!  there  are  an  upper  and  a 
lower  values  of  Lm  P  (joy )  vs  to,  and  ZP,  vs  co, ,  respectively. 


Figure  4:  Bode  Plots  of  6  LTI  Plants:  the  Range  of  Parameter  Uncertainty. 

Control  System  Performance  Specifications  (P.S.)  —  The  desired  system  output,  y(t),  is  to  lie  between  the 
specified  upper  and  lower  bounds  y(t)u  and  y(t)L,  respectively,  as  shown  in  Fig.  5(a).  The  performance 
figures  of  merit  (FOM),  based  upon  a  step  input  signal  r(t)  =  R0  u./(t),  are  indicated  in  Fig.  5(a).  The  FOM 
are:  Mp  (peak  overshoot),  tr  (rise  time),  tp  (peak  time),  and  ts  (settling  time).  The  corresponding  P.S.  in  the 
frequency  domain  are  By  and  BL,  the  upper  and  lower  bounds  respectively.  The  P.S.  are:  peak  overshoot 
(Lm  Mm)  and  frequency  bandwidth  (oy, )  as  indicated  in  Fig.  5(b).  [Note:  the  increasing  Or  (joy)  above  the  0  dB 
crossing  characteristic.]  With  negligible  sensor  noise,  sufficient  control  effort  authority,  and  a  stable 
linear-time-invariant  {LTI)  minimum-phase  ( m.p .)  plant  a  LTI  compensator  may  be  designed  to  achieve 
the  desired  control  system  P.S. 
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t 


(a) 


(b) 


Figure  5:  Desired  System  Performance  Specifications:  (a)  Time  Domain  Response  Specifications; 
(b)  Frequency  Domain  Response  Specifications. 


QFT  Design  Overview  —  The  QFT  design  objectives  are  achieved  by  representing  the  characteristics  of 
the  plant  and  the  desired  system  P.S.  in  the  frequency  domain.  These  representations  are  used  to  design  a 
compensator  (controller)  that  results  in  achieving  the  desired  P.S.  The  nonlinear  plant  characteristics  are 
represented  by  a  set  of  J  LTI  transfer  functions  that  cover  the  range  of  structured  parametric  uncertainty. 
Two  LTI  transfer  functions  that  yield  the  system’s  desired  P.S.  are  synthesized  which  result  in  the  upper 
By  and  lower  BL  boundaries  for  the  design  as  shown  in  Fig.  5.  The  effect  of  parameter  uncertainty  is 
reduced  by  shaping  the  open-loop  frequency  responses  so  that  the  Bode  plots  of  the  J  closed-loop  systems 
fall  between  the  boundaries  of  B(;  and  BL,  while  simultaneously  satisfying  all  P.S.  The  normal  open- loop 
transfer  function  L0  is  shaped  so  that  the  stability,  tracking,  disturbance,  and  cross-coupling  (for  MIMO 
systems)  boundaries  (bounds)  on  the  Nichols  chart  (N.C.)  are  satisfied  which  in-tum  results  in  achieving 
the  desired  P.S. 

QFT  Basics  —  Consider  the  Control  system  of  Fig.  6  where  G(s)  is  a  compensator,  F(s)  is  a  prefilter,  and 
IP  is  a  nonlinear  plant.  To  accomplish  a  QFT  design  for  this  system  the  nonlinear  plant  is  described  by  a  set 
of  J  m.p.  LTI  plants,  i.e.,  IP  =  (P,(s)}  ( i  =  1,2,...,J),  which  define  the  structured  plant  parameter 
uncertainty.  Note  that  for  MIMO  systems  the  elements  py  of  the  mxm  plant  matrix  P  can  be  n.m.p. 
For  MISO  systems  the  discussion  is  restricted  to  m.p.  plants.  The  magnitude  variation,  8P(je0i),  is  due  to 
the  plant  parameter  uncertainty  which  is  depicted  by  the  Bode  plots  of  the  LTI  plants  (see  Fig.  7)  for  the 
example  of  Figs.  3  and  4.  The  J  data  points  [log  magnitude  (Lm)  and  phase  angle  (<[>)],  for  each  value  of 
frequency,  co  =  C0i,  are  plotted  on  the  N.C.  The  contour  that  is  drawn  through  the  data  points,  for  co;, 
describes  the  boundary  of  the  region  that  contains  all  J  points.  The  contours  are  referred  to  as  templates 
which  represent  the  region  of  structured  plant  parametric  uncertainty  on  the  N.C.  and  are  obtained  for 
specified  values  of  frequency,  co  =  C0i,  within  the  bandwidth  (BW)  of  concern.  The  six  data  points 
(Lm  and  (|))  for  the  example  of  Fig.  3,  for  each  value  of  cox,  are  obtained  as  shown  in  Fig.  8a.  These  points 
are  used  to  plot  the  templates,  for  each  value  of  co,  as  shown  in  Fig.  8b  for  co,  =  3  rad/sec  and  in  Fig.  8c  for 
other  values  of  co;.  The  system  P.S.  are  represented  by  LTI  transfer  functions  whose  corresponding 
Bode  plots  are  shown  in  Fig.  7  by  the  upper  and  lower  bounds  By  and  BL,  respectively. 


Figure  6:  Compensated  Nonlinear  System. 
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Figure:  7  LTI  Plants. 
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Figure  8:  (a)  Bode  Plots  of  6  LTI  Plants;  (b)  Template  Construction  for  ©i  =  3  rad/sec; 
(c)  Construction  of  the  Nichols  Chart  Plant  Templates. 


1  -6 


RTO-EN-SCI-142 


Quantitative  Feedback  Technique  (QFT):  Bridging  the  Gap 


QFT  Design  —  The  tracking  design  objective  is  to  synthesize  a  compensator  G(s)  of  Fig.  6  which  results 
in  satisfying  the  desired  P.S.  of  Fig.  5  and  in  the  corresponding  closed-loop  frequency  responses 
TL(  shown  in  Fig.  9.  This  results  in  the  SL(ja>i)  of  Fig.  9,  of  the  compensated  system,  being  equal  to 

or  smaller  than  8P(j00i)  of  Fig.  7  for  the  uncompensated  system  and  that  it  is  equal  or  less  than  8R(jcOi), 
for  each  value  of  CO;  of  interest;  that  is: 


§L(j»i)  <  5r  (Jcoi)  <5P(jc0i) 


The  next  step  is  to  synthesize  a  prefilter  F(s)  of  Fig.  6  which  results  in  shifting  and  reshaping  the  Tl  ; 
responses  in  order  that  they  lie  within  the  By  and  BL  boundaries  in  Fig.  9  as  shown  in  Fig.  10. 


Figure  9:  Closed-loop  Responses: 
LTI  Plants  with  G(s). 


FigurelO:  Closed-loop  Responses: 
LTI  Plants  with  G(s)  and  F(s). 


Therefore,  the  QFT  robust  design  technique  assures  that  the  desired  performance  specifications  are 
satisfied  over  the  prescribed  region  of  structured  plant  parametric  uncertainty. 


Insight  to  the  QFT  Technique 

Fig.  6  is 


The  open-loop  plant  transfer  function  for  the  position  control  system  of 


P,(s)=  Ka  =  K' 
s(s  +  a)  s(s  +  a) 


(7) 


where  K'  =  Ka  and  i  =  1,2,...,  J.  The  Lm  variation  due  to  the  plant  parameter  uncertainty,  for  J  =  6, 
is  depicted  by  the  Bode  plots  in  Fig.  4.  The  loop  transmission  L(s)  of  the  unity-feedback  system  of  Fig.  6 
is  defined  as 

L,(s)  =  G(.v)P,(.s)  (8) 


The  overall  system  control  ratio  TR,  which  includes  the  prefilter,  is  given  by: 


Tr,(s)  = 


F(s)L,(s) 
1  +  L,  (s) 


(9) 


for  each  i  plant  and  whose  control  ratio  TL,  is: 


Y  L, 
Rl  1  +  L, 


(10) 
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Results  of  Applying  the  QFT  Design  Technique  —  Proper  application  of  the  robust  QFT  design 
technique  requires  the  utilization  of  the  prescribed  P.S.  from  the  onset  of  the  design  process,  the  selection 
of  a  nominal  plant  P0  from  the  J  LTI  plants,  and  the  proper  loop  shaping  of  L0(s)  =  G(s)P0(s)  which  results 
in  obtaining  a  synthesized  G(s)  that  satisfies  the  desired  P.S.  The  last  step  of  this  design  process  is  the 
synthesis  of  a  prefilter  that  ensures  that  the  Bode  plots  of  TR;  all  lie  between  the  upper  and  lower  bounds 
By  and  BL. 


Insight  to  the  Use  of  the  Nichols  Chart  (N.C.)  in  the  QFT  Technique  —  A  review  of  the  use  of  the 

N.C.  and  an  insight  as  to  how  it  applies  to  the  QFT  technique  is  now  provided. 

Open-Loop  Characteristics  -  For  the  nominal  plant  P0(jco),  the  nominal  loop  transmission  function  is 

Lm  L0  =  Lm  GP0  =  Lm  G  +  Lm  P0  (11) 

Whereas  for  all  other  plants,  P(joj),  the  loop  transmission  function  is 

Lm  L  =  Lm  GP  =  Lm  G  +  Lm  P  (12) 


Thus,  for  co  =  cOj,  the  variation  SP(j(Oi)  in  Lm  L(j oji)  is  given  by 

Sp(ja>i)  =  Lm  L(jcOi)  -  Lm  L0(jec>i)  =  Lm  P(jcOi)  -  Lm  Po(jc0i)  (13) 

and  its  phase  angle  variation  is  given  by 

ZAP(jco.)  =  ZL  -  ZL0  =  (ZG  +  ZP)  -  (ZG  +  ZP0)  =  ZP  -  ZP0  (14) 

The  expression  Lm  P(]ec>i)  =  Lm  Po(jc0i)  +  SP(jcOi),  obtained  from  Eq.  (13),  is  substituted  into  Eq.  (12)  to 
yield 

Lm  L(jcOi)  =  LmG(jco;)  +  LmPo(je0i)  +  SP(jcc>i)  (15) 


where  SP(j(Oi),  is  the  variation  from  the  nominal  point. 


Closed-Loop  Characteristics  -  The  closed-loop  system  characteristics  are  obtained,  for  a  given  G(jco)  and 
P0(jco),  from  the  plot  of  Lm  L0(jco)  vs.  ZL0  shown  on  the  N.C.  in  Fig.  11.  Also,  is  shown  a  plot  of  a  template, 
TTP(j<x>i),  whose  contour  is  based  upon  the  data  obtained  for  co  =  to,  from  Fig.  4.  The  template  represents  a 
region  of  plant  parameter  uncertainty  for  <x>i  as  expressed  mathematically  by  Eqs.  (14)  and  (15).  From  the 
loop  transmission  plot  and  its  intersections  with  the  M-  and  a-contours  (note:  the  a  contours  are  not 
shown),  the  closed-loop  frequency  response  data  may  be  obtained  for  plotting  M0  and  a  vs.  co.  In  Fig.  12 
a  plot  of  M0  vs.  co  is  shown,  where 


M0(jcy)Z«(jcy) 


R(jco) 


L0(jft>) 

1  +  L0(j&>) 
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Figure  11:  Nominal  Loop  Figure  12:  Closed-loop  Responses  obtained  from 

Transmission  Plot  with  Fig.  11.  Plant  Parameter  Area  of  Uncertainty. 

Parametric  Variation  N.C.  Characteristics  -  As  an  example  consider  that  for  Eq.  (8)  G  =  1  Zl)°.  If  point 
A  on  the  template  in  Fig.  1 1  represents  Lm  P0  vs.  AP0,  then  a  variation  in  P  results  in  a  horizontal 
translation  in  the  angle  of  P,  given  by  Eq.  (14)  and  in  a  vertical  translation  in  the  Lm  value  of  P,  given  by 
Eq.  (15).  The  translations  are  shown  in  Fig.  11  at  points  B,  C,  and  D.  The  variation  8P(j©i)  of  the  plant, 
and  in  turn  L(jcOi),  from  the  nominal  value  for  co  =  w,,  over  the  range  of  plant  parameter  variation, 
is  described  by  the  template  JPO©!)  shown  in  Fig.  11.  For  example,  point  A  on  the  template  represents  the 
nominal  plant  P  =  Po(jc0i).  Its  corresponding  closed-loop  frequency  response  is: 

Lm  Ma  =  Lm  P  =  -6  dB  (16) 

For  P(jco)  =  P,(j©):  point  C  in  Fig.  1 1,  its  corresponding  closed-loop  frequency  response  is 

Lm  Mc  =  Lm  a  =  -2  dB  (17) 

These  values  are  plotted  in  Fig.  12.  Note  that  point  A  represents  the  minimum  value  of  Lm  M(j©;)  at 
co  =  CO;,  i.e., 

Lm  Tmin  =  [Lm  M(jcOi)]min  =  Lm  P  (18) 


Point  C  represents  the  largest  value  of  Lm  M(jo),)  at  ©  =  ©,,  i.e., 

Lm  Tmax  =  [Lm  M(j©i)]max  =  Lm  a  (19) 

For  the  range  of  plant  parameter  variation  described  by  the  template  the  maximum  variation  in  Lm  M, 
denoted  by  SL(j©i),  for  this  example  is 

SL  G  )  -  Lm«  -  Lm  J3  =  -  2  -  (-6)  =  4  dB 


When  L0(s)  is  properly  synthesized  then  8[.(j©j)  <  SR(j©i).  Shown  in  Fig.  12  is  point  E,  midway  between 
points  A  and  C,  that  corresponds  to  point  E  in  Fig.  1 1  which  lies  within  the  variation  template  fiP(j©j). 
This  procedure  is  repeated  to  obtain  the  maximum  variation  5L(j©i)  for  each  value  of  frequency  ©  =  w, 
within  the  desired  BW  (see  Fig.  12).  From  this  figure  it  is  possible  to  determine  the  variation  in  the  control 
system’s  figures  of  merit  due  to  the  plant’s  parameter  uncertainty.  This  graphical  description  of  the  effect 
of  plant  parameter  uncertainty  on  the  system’s  performance  is  the  basis  of  the  QFT  technique. 
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Benefits  of  QFT  —  A  design  is  achieved  that  is  robust  for  the  full  operating  envelope  and  is  insensitive  to 
structured  plant  parameter  variation.  QFT  also  has  the  following  features:  the  design  limitations  are 
apparent  up  front  and  during  the  design  process,  the  achievable  P.S.  can  be  determined  early  in  the  design 
process,  the  designer  is  able  to  quickly  redesign  for  the  changes  in  P.S.  with  the  aid  of  a  QFT  CAD 
package,  and  the  structure  of  the  compensator  (controller)  is  determined  up  front.  Thus,  these  benefits 
result  in  less  development  time  for  a  full  envelope  design. 

Summary  —  This  section  presents  the  concepts  of  the  QFT  technique  in  such  a  manner  that  one  can 
readily  grasp  the  fundamentals  and  appreciate  its  transparency  in  bridging  the  gap  between  theory  and  the 
real  world.  For  a  more  theoretical  discussion  of  QFT,  such  as  existence  theorems  and  nonminimum-phase 
(n.m.p.)  plants,  the  reader  is  referred  to  the  literature  listed  in  the  reference  section  of  Floupis,  et  al., 
(1999).  It  is  highly  recommended  that  the  reader  read  the  excellent  paper  by  Horowitz,  (1991). 


1.2  MISO  Analog  Control  System 


Introduction  —  As  an  introduction  to  the  MISO  design  first  consider  the  mxm  control  system  of  Fig.  13 
which  is  shown  later  can  be  represented  by  the  equivalent  rrf  MISO  systems  of  Fig.  14.  This  MISO 
structure  representation  of  a  MIMO  system  forms  the  basis  for  the  development  of  the  QFT  design 
technique. 


Figure  13:  MIMO  Feedback  Structure. 


Figure  14:  Effective  MISO  Loops  Two-by-Two 
(boxed  in  loops)  and  Three-by-Three  (all  9  loops). 


MISO  Control  System  —  The  MISO  QFT  design  technique  is  presented  for  control  systems  having 
minimum-phase  (m.p.)  LTI  plants.  The  control  ratios  for  tracking  (D  =  0)  and  disturbance  rejection 
(R  =  0)  for  the  system  of  Fig.  15  are,  respectively: 

T  N  F(s)G(s)P(s)  F(s)L(s) 

Tr  (s)  = - = -  (20) 

1  +  G(s)P(s)  1  +  L(s) 


Td(s)  = 


P(s) 

1  +  G(s)P(s) 


P(s) 

1  +  L(s) 


(21) 


R^O- 


D 

A 


O- 


9 


Y 


Figure  15:  Compensated  MISO  Control  Systems. 
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The  robust  design  objective  is  to  design  G(s)  and  F(s)  of  Fig.  15  for  the  plant  having  parameter 
uncertainty.  The  design  procedure  to  achieve  this  objective  is  as  follows: 

Step  1:  Synthesize  tracking  models 

Step  2:  Synthesize  disturbance  model 

Step  3 :  Specify  J  LTI  plant  models  (define  uncertainty) 

Step  4:  Obtain  templates  at  specified  eo,  (describes  uncertainty) 

Step  5 :  Select  nominal  plant  P0(s) 

Step  6:  Determine  stability  contour  (U-contour)  on  N.C. 

Steps  7-9:  Determine  disturbance,  tracking,  and  optimal  bounds 

Step  10:  Synthesize  nominal  L0(s)  =  G(s)P0(s)  to  satisfy  all  bounds  and  the  stability  contour  in  order  to 

obtain 

G(s)  =  Lo(s)/P0(s) 


Step  1 1 :  Synthesize  prefilter  F(s). 

Step  12:  Simulate  linear  system  (J  time  responses) 

Step  13:  Simulate  with  nonlinearities 


Synthesize  Tracking  Models  —  Based  upon  the  FOM  specifications  [see  Fig.5a]  MP,  tp,  ts,  and  tr,  for  a 
unit  step  forcing  function,  underdamped  and  overdamped  tracking  models  are  synthesized. 

Simple  Second-order  Models  -  Models  are  synthesized  for  an  upper  bound  y(t)i:  and  for  a  lower 
bound  y(t)L  that  yield  bounds  between  which  the  acceptable  response  y(t)  must  lie.  For  m.p.  plants, 
only  tolerance  on  |TR(jec>j)|  need  be  satisfied  where  for  nonminimum-phase  (n.m.p.)  plants  the  tolerances 
on  ZfTjdjeOi)]  must  also  be  satisfied. 

Synthesized  tracking  control  ratios  -  The  tracking  models  are: 

(ft>n/a)(s  +  a) 

TrU  (s)  -  s2  +  2gw„S  +  W„  UPPerB°“"d  <22> 


Trl(s) 


K 

Lower  Bound 

(s  -  cr i )(s  -  cr2)(s  -  cr3) 


(23) 


These  models  are  synthesized  so  that  in  Fig.  5b  for  co;  >  eoh  SR(j  00i)  increases  in  order  to  simplify 
the  process  of  synthesizing  L0(s)  =  G(s)P0(s).  Note  that  the  frequency  bandwidth  (BW),  0  <  co  <  (Oh 
[see  Fig.  5b]  is  determined  by  the  -12  dB  line. 

Disturbance  Models  —  Based  upon  a  unit  step  disturbance  input  (D  of  Fig.  15)  the  simplest  disturbance 
model  is: 

|TD(jcq)|  =  |Y(jco)/D(jco)|  <  ap  (maximum  allowable  magnitude  of  |TD|) 
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Thus,  the  frequency  domain  specification  is: 

Lm  TD(jeo)  <  Lm  ap  over  the  BW 

This  results  in  only  an  upper  bound. 

J  LTI  Plant  Models  —  The  following  simple  plant  is  chosen  for  illustrative  purposes: 

Ka  K' 

p  rs\  = _ = _ 

s(s  +  a)  s(s  +  a) 


(24) 


where  K'  =  Ka,  K  e  {1,  10}  and  a  e  { 1 ,  10}.  The  region  of  plant  parameter  uncertainty  (see  Fig.  16)  which 
is  described  by  the  J  LTI  plants  defined  by  points  1,  2,  3,  4,  5,  and  6.  These  models  define  the  region  of 
plant  parameter  uncertainty  as  shown  in  the  figure. 


K 


2  3 

■t - • 


Region  of  plant 
parameter  uncertainty 


T  + 


5 


0 


a 


Figure  16:  Region  of  Plant  Parameter  Uncertainty. 


Plant  Templates  of  P„  3P(j(0j)  —  The  overall  tracking  control  ratio  of  Fig.  15,  where  L  =  GP,  is 


Tr  =F 


L 

1  +  L 


which  yields 


LmTR 


=  Lm  F  +  Lm 


L 

1  +  L 


Subtracting  Lm  F  from  sides  of  the  above  equation  yields: 


A(LmTR)  =  LmTR  -  Lm 


L 

1  +  L 


The  change  in  Tr  due  to  the  uncertainty  in  P  (F  is  LTI),  is: 


Sr  =  A(LmTR( 


Lm  F)  =  Lm 


L 

1  +  L 


(25) 


for  i  =  1,2,  ...  ,  J  plants.  The  nominal  loop  transmission  L0  and  the  prefilter  F  are  designed  (synthesized) 
so  that  the  actual  value  of  Lm  T  always  lies  between  By  and  BL  of  Fig.  5b.  Thus,  draw  the  templates  that 
characterize  the  variation  of  plant  uncertainty  onto  the  Nichols  Chart  (N.C.)  for  various  values  of  oj,  over 
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the  BW.  The  six  points  on  the  boundary  of  the  region  of  plant  uncertainty  are  mapped  onto  the  N.C.  in 
Fig.  17  where  1  — »  A,  2  — »  B,  4  — >  C,  5  — »  D.  Note  that  point  3  lies  between  B  and  C  and  point  6  lies 
between  D  and  A.  A  curve  is  drawn  through  these  points,  as  indicated  in  the  figure,  and  the  shaded  area  is 
labeled  3P(j  1).  Templates  for  other  values  of  co,  are  obtained  in  the  same  manner. 


Figure  17:  N.C.  Characterizing  Eq.  (7)  over  the  Region  of  Uncertainty. 

Characteristic  of  Templates  —  The  templates  of  Fig.  8b  have  the  following  characteristics:  (1)  starting 
from  low  values  of  co;  the  angular  width  of  the  templates  become  larger;  (2)  for  increasing  values  of 
©i  the  templates  then  become  narrower;  and  (3)  they  eventually  approach  a  straight  line  of  height  V  dB 
[see  Eq.  (27)]. 

Nominal  Plant  —  Any  of  the  J  plants  can  be  chosen  as  the  nominal  plant  PG.  One  can  select, 
if  specifications  permit,  a  plant  whose  N.C.  point  is  always  at  the  lower  left  comer  for  all  templates. 
In  Houpis,  et  al.  ( 1 999)  general  guidelines  are  given  in  Chapter  9  for  selecting  the  nominal  plant. 

U-Contour  (Stability  Bound)  —  The  value  of  Mm  (see  Fig.  5b)  identifies  the  desired  maximum 
overshoot  for  a  closed-loop  system.  Thus,  from  this  figure  the  bound  MP  =  Mm  is  chosen.  This  bound 
identifies  the  M-contour  on  the  N.C.  Mm  =  ML  (see  Fig.  5).  This  contour  establishes  a  region  that  can 
not  be  penetrated  by  the  templates  and  by  L0(jco).the  dominating  constraints  on  L0(jco).  The  top  portion 
(efa)  of  the  ML  contour  becomes  part  of  the  U-contour  (See  Fig.  18).  The  limiting  value  of  the  plant 
transfer  function,  for  a  large  problem  class,  approaches 

lim[P(j®)]=—  (26) 

CO 

where  X  is  the  excess  of  poles  over  zeros  of  P(s).  The  plant  template,  for  this  problem  class,  approaches  a 
vertical  line  of  length  equal  to 

A  =  lim  [LmPmax  -  LrnPmin]  =  LrnKmax  -  LmKmin  =  V  dB  (27) 

6)—>cC 

The  nominal  plant,  for  this  example,  is  chosen  at  K  =  Kmin  (point  1  of  Fig.  16)  which  results  in  the 
U-contour  abcdefa  of  Fig.  18. 
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Figure  18:  U-contour  Construction  (Stability  Contour). 

Optimal  (Composite)  Bounds  B0(jcOi)  on  L0(ja>j)  —  By  applying  the  procedures  given  by  Houpis,  et  al., 
(1999),  the  tracking  bounds  BR(jec>i),  and  the  disturbance  bounds  BD(jec>i)  are  determined  which  in  turn  result 
in  obtaining  the  composite  bounds  B0(jo0i)  that  become  the  bounds  required  in  synthesizing  L0(  jco;). 

Tracking  Bounds  on  La  -  The  determination  of  the  tracking  bounds  BR(jcOi)  requires  that 

SLOi)  =  ATrOO  <  5R(jcoi)  dB 

be  satisfied  for  all  L,(ja>i)  [see  Figs.  5b,  9  and  10].  The  procedure  for  determining  these  bounds  is  based  on 
a  nominal  plant  L0,  is  as  follows:  at  each  specified  value  eo;  and  the  use  of  the  corresponding  template 
(or  CAD),  see  Fig.  19,  the  nominal  point  A  on  the  template  is  moved  up  or  down  along  the  N.C.  angle 
grid  line  without  rotating  the  template;  this  template  is  moved  until  it  is  tangent  to  two  M-contours  whose 
difference  in  their  M  values  is  essentially  equal  to  SR.  When  this  condition  is  achieved,  the  location  of  the 
nominal  point  on  the  template  becomes  a  point  on  the  bound  BR(jec>i).  This,  procedure  is  repeated  on 
sufficient  angle  grid  lines  to  provide  enough  points  to  draw  BR(jeo,)  and  is  repeated  for  all  values  of 
frequency  for  which  templates  have  been  obtained. 


f. 


Figure  19:  Graphical  Determination  of  BrGco,). 
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Disturbance  Bounds  Procedure  (MISO  System)  -  Two  cases  of  disturbance  bounds  are 
considered  as  indicated  in  Fig.  20.  These  cases  are:  Case  1  [d2(t)  =  D0u_i,(t),  di(t)  =  0]  and  Case  2 
[di(t)  =  D0u_i,(t),  d2(t)  =  0]. 


Figure  20:  A  Feedback  Structure. 


Case  1:  From  Fig.  20,  the  disturbance  control  ratio  for  input  d2(t)  is 

1 

To  (s)  = -  (28) 

1  +  L 


Substituting  L  =  1  /(  into  Eq.  (28)  yields: 

Td  (s)  =  ~  (29) 

1  +  f 

that  has  the  mathematical  format  required  to  use  the  N.C.  Over  the  specified  BW  it  is  desired  that 
|TdG®)I  «  1,  which  results  in  the  requirement,  from  Eq.  (28),  that  |L(jco)|  »  1  (or  |f(j©)|  «  1),  i.e., 

|TD0)|*i — - — r  =  KH  (30) 

|L(j<a?)| 


A  time-domain  tracking  response  based  upon  r(t)  =  u.i(t)  often  specifies  a  maximum  allowable  peak 
overshoot  Mp,  which  in  the  frequency  domain  this  specification  may  be  approximated  by: 


Mr  (jo)  =  TR(jco) 


YQco) 

R(jco) 


<  Mm  «  Mp 


(31) 


The  corresponding  time-  and  frequency-domain  disturbance  response  characteristics,  based  upon  the  step 
disturbance  forcing  function  d2(t)  =  u_i(t),  are  determined  in  the  same  manner,  i.e., 


mD(t) 


y(t) 

d(t) 


<  ap 


for  t  >  tx 


(32) 


and 


|md  (jco)|  =  |td  ( jco) 


Y(j©) 

D(jra) 


<  am 


ap 


(33) 


The  disturbance  control  ratio  of  Eq.  (29)  must  have  the  desired  characteristic  over  the  desired  BW  0  <  ©  <  ©2 
for  which  |  L(j©)  |  »  1  and  |  f(jeo)  |  «  1.  Thus,  Eq.  (31)  applies  within  this  BW  region.  Consider  for  a 
given  plant  L(s)  that  its  plot  of  Lm  L  vs  ZL  is  shown  in  Fig.  21.  The  plot  is  tangent  to  the  Lm  M  =  0.5  dB 
contour  at  ©  =  ©m.  Table  1  contains  data  for  two  points  on  the  Nichols  plot  of  Fig.  21.  Also,  shown  is  the 
plot  of  Lm  f(j©)  vs.  Zb  for  these  two  points. 
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Table  1 :  Data  Points  for  a  Nichols  Plot 


CO 

Lm  L 

ZL 

Lm  £ 

Z£ 

CDi 

21 

-96° 

-21 

96°  (or  -264°) 

C02 

15 

O 

OO 

o\ 

i 

-15 

98°  (or  -262°) 

Consider  that  the  template  AP(ja>i)  for  a  given  plant  P  having  uncertain  parameters,  for  a  given  o)„  has 
equal  dB  differences  along  its  A-B  and  C-D  boundaries,  i.e.,  for  a  given  ZP(jec>j), 

A(LmPB  -  LmPA)  =  A(LmPc  -  LmPD)  =  10  dB 

This  template  is  arbitrarily  set  on  the  N.C.  as  shown  in  Fig.  21.  The  data  corresponding  to  the  template 
location  shown  in  Fig.  21  is  given  in  Table  2.  The  plot  of  Lm  £(jeo)  vs.  Zl  for  these  2  points  is  also 
shown  in  this  figure.  The  template  of  the  reciprocal,  Lm  [1/P(jec>i)],  is  arbitrarily  set  on  the  NC  in  Fig.  21 
for  the  same  frequency  as  for  the  template  of  Lm  P(joOi)  and  for  the  angles  of  Table  1.  Note,  the  template 
of  Lm  [1/P(]ec>i)]  is  the  same  as  the  template  of  Lm  P(jtOi)  but  is  rotated  by  180°.  It  is  located  by  first 
reflecting  the  template  of  Lm  P(jec>i)  about  the  -180°  axis,  “flipping  it  over”  vertically,  and  then  moving 
it  up  or  down  so  that  it  lies  between  -5  and  -20  dB.  For  the  arbitrary  location  of  the  template  of 
Lm  [1/P(jffli)],  note  that: 

(a)  Pab  =  180°  +  ZP  =  180°  -  120°  =  60° 

pCD  =  180°  +  ZP  =  180°  -  80°  =  100° 

where  for  1/P,  the  corresponding  angles  are 
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(b)  The  templates  of  Lm  P(j co, )  are  used  for  the  tracker  case:  T  =  L/(l  +  L)  where  L  =  GP. 

(c)  The  templates  of  Lm  [1/P(]ec>i)]  used  for  the  disturbance  rejection  case:  TD  =  1/(1  +  L)  =  £/(  1  +  £) 
where  £  =  1/GP  =  1/L. 


Table  2:  Data  Points  for  the  Templates  of  Fig.  21 


Points 

ZP 

Zl/P 

A,B 

-120° 

120°  (or -240°) 

C,D 

-80° 

80°  (or  -280°) 

Since  Lm  L  is  desired  for  the  disturbance-rejection  case  of  Eq.  (28)  then  Eq.  (29)  must  be  utilized  to 
determine  the  corresponding  disturbance  bounds.  Thus,  the  N.C.  of  Fig.  21  is  rotated  180°  and  is  redrawn  in 
Fig.  22  in  order  to  determine  the  disturbance  boundary  BD(ja>j)  for  L(jcOi)  =  l/£{j(oi)  on  this 
rotated  N.C.  Since  Lm  £(jco)  =  Lm  [ l/L(joo)]  =  -  Lm  L(jeo),  then  a  negative  value  of  Lm  £  yields  a 
positive  value  for  Lm  L  as  shown  in  Fig.  22.  Since  L  =  KL'  =  Vi  ,  then 

£(]co)  -  K_1f  (ja>)  (34) 


(Note  Use  the  negative  angle  for/tL  since  n  >  w) 


Lm  f  Lm  L 

-24-p24 
-20- -20 
-16--16 
-12--12 
-8  --  8 
-4-4  3 

o  -  o  ! 

c 

c 

4  --4  | 
8  ---8 
12  -  —  12 
16  -  —  16 
20  -—20 


24 


-24 


28  -1-28 


Figure  22:  Rotated  Nichols  Chart. 


If  r(jco)  is  given  and  it  is  required  to  determine  K'1  to  satisfy  Eq.  (33),  then  the  plot 
Lm  f'(joo)  vs.  Zi'( jeo)  must  be  raised  or  lowered  until  it  is  tangent  to  the  Lm  ap-contour  (|TD|max  =  am). 
The  amount  A  by  which  the  plot  is  raised  or  lowered  yields  the  value  of  K,  i.e.,  Lm  K’1  =  A.  Note,  that  this  is 
the  same  procedure  used  for  obtaining  the  tracking  bounds,  except  that  the  adjustment  in  Lm  f(jeo)  is  K"1. 

The  rotated  N.C.  is  used  to  determine  directly  the  disturbance  boundaries  BD(jcOi)  for  LD(jc0i).  Point  A  for 
the  simple  plant  of  this  design  example  corresponds  to  the  nominal  plant  parameters  and  is  the  lowest 
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point  of  the  template  ffP(jcOi).  This  point  is  again  used  to  determine  the  disturbance  bounds  BD(jo)i). 
The  lowest  point  of  the  template  must  be  used  to  determine  the  bounds  and,  in  general,  may  or  may  not  be 
the  point  corresponding  to  the  nominal  plant  parameters.  Based  upon  Eqs.  (28)  -  (30) 

-  LmTo  =  Lm[l  +  L]  >  -  Lm«(jmi)  >  0  dB  (35) 

where  a(ja>i)  <  0.  Since  |L|  »  1  in  the  BW,  then 

-LmTo  =  LmL  >  -Lma(jcoi)  =  -  JdOcoO  (36) 


In  terms  of  L(jeOj),  the  constant  M-contours  of  the  N.C.  can  be  used  to  obtain  the  disturbance  performance 
Td.  This  requires  the  change  of  sign  of  the  vertical  axis  in  dB  and  the  M-contours,  as  shown  in  Fig.  22. 
The  procedure  for  determining  the  boundaries  BD(jeOi)  is  given  by  Houpis,  et  al. ,  (1999). 


Case  2  [di(t)  =  D0u_i(t),  d2(t)  =  0]:  The  disturbance  control  ratio  for  this  case  (see  Fig.  20)  is: 


Td  (j  g>) 


1  +  G(ja>)P(ja>) 


(37) 


This  equation,  assuming  point  A  of  the  template  (Fig.  23)  represents  the  nominal  plant  P0,  is  multiplied  by 
PG/P0  and  is  rearranged  to  obtain: 


Td 


1 


1 

—  +  G 

P 


P  P 

~  +  GP0  —  +  Lc 

P  P 


w 


where  W  =  (Pq/P)  +  L0.  From  Eq.  (38),  set  Fm  TD  =  SD  =  Fm  ap  and  obtain 

Fm  W  =  Fm  P0  -  SD 


(38) 


(39) 


Templates  are  drawn  in  polar  coordinates  for  each  value  of  a>i.  Shown  in  Fig.  23  is  a  template  for  eo  =  cc>i. 
From  Eq.  (39)  evaluate  IWOca)!  for  eo  =  00;.  This  magnitude  is  used  in  conjunction  with  the  following 
equation 

W(jcoi)  =  [P0(jcoi)/P(jcoi)]  +  L0(j®i) 

to  obtain  a  graphical  solution  for  Bq(jeoj)  (see  Fig.  24).  Thus,  the  N.C.  bounds  for  the  specified  values  of  CO; 
are  given  by: 

BD(jo>i)  =  Lm  Bd(jeOi)  and  Z  BD(ja>i)  =  -  180° 
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Figure  24:  Graphical  Evaluation  of  Bd(jcoi). 


Optimal  Bounds  -  The  optimal  bound  B0(jo>i)  for  the  case  shown  in  Fig.  25a  is  composed  of  the  portions 
of  brg  <x>i)  and  BD(jcOi)  that  have  the  largest  dB  values.  For  the  case  shown  in  Fig.  25b  the  outermost 
portions  of  BR(jec>j)  and  BD(jeOj)  becomes  the  perimeter  of  B0(j(Oi)  (yields  the  greatest  interior  area). 
For  example,  the  synthesized  L0(ju)i)  must  lie  on  or  just  above  the  bound  B0(jcOi)  of  Fig.  25  a. 


Figure  25:  Composite  B0(j»i). 


Synthesizing  (or  Loop  Shaping)  —  The  shaping  of  L0(jco),  the  dashed  curve  shown  in  Fig.  26  is  achieved 
as  follows:  (1)  a  point  such  as  Lm  L0(]2)  must  be  on  or  above  the  U-contour  and  BR(j2);  (2)  to  satisfy  the 
performance  specifications  L0(jec>i)  cannot  violate  (intersect)  the  U-contour;  (3)  L0  should  closely  follow 
the  U-contour  up  to  40  r/s  and  must  stay  below  it  above  40  r/s;  and  for  this  example  it  must  be  at  least  a 
Type-1  function  with  at  least  one  pole  at  the  origin.  In  synthesizing  the  rational  function  L()(.v)  that  satisfies 
the  above  involves  building  up  the  function: 


W 

LoG®)=  Lok(jfi>)  =  Po (jco)  n[KkGk(j®)] 


k=0 


(40) 


W 

where  for  k  =  0,  G0  =  I  Z0°,  and  K  =  n  Kk  .  In  order  to  minimize  the  degree  of  G(s)  assume  initially  that 

k  =  0 

L0(jeo)  =  P0(]eo)  in  Eq.  (40)  and  then  build  up  L0(jeo)  term-by-term  so  that  the  point  L0(jcOi)  lies  on  or 
above  the  corresponding  B0(ja>i)  and  it  stays  just  outside  the  U-contour.  The  design  of  a  proper  L0(s) 
guarantees  only  that  the  variation  SL  is  less  than  or  equal  to  SR(ja>i).  Once  a  proper  L0(s)  is  synthesized  then  the 
compensator  is  given  by  G(s)  =  Lo(s)/P0(s).  For  this  example  L0(jeo)  slightly  intersects  the  U-contour.  Since 
the  QFT  technique  has  the  inherent  characteristic  of  resulting  in  an  over-design,  thus  for  the  first  trial 
design  no  effort  is  made  to  fine  tune  G(s).  If  the  simulation  results  are  not  satisfactory  then  L0  is  fine 
tuned. 


RTO-EN-SCI-142 


1  - 19 


Quantitative  Feedback  Technique  (QFT):  Bridging  the  Gap 


ORGANIZATION 


/ 


Figure  26:  Shaping  of  L0(jco)  on  the  Nichols  Chart  for  the  Plant  of  Eq.  (23). 


Prefilter  Design  F(s)  —  The  design  of  L0(s)  guarantees  only  that  the  variation  in  |TR(]eo)|  is: 

LmTR(j(o)max,  -  Lm  TR(j(o)min  <  5R(jeo) 


Thus,  the  design  procedure  involves  positioning 


Tl,  (s)  = 


L/(s) 


1  +  L,(s) 

so  that  it  lies  between  Bu  and  BL  (see  Fig.  5b)  for  all  J  plants.  The  filter  design  procedure  is  as  follows: 


(41) 


Step  1:  Place  the  nominal  point,  for  this  example  Point  A,  of  co;  on  the  plant  template  on  to  the 

L0(j©i)  point  of  the  L0(ja>)  curve  (see  Fig.  27). 


Figure  27:  Prefilter  Determination. 
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Step  2: 

Step  3 : 


Traversing  the  template,  obtain  the  values  of  Lm  Tmax  and  Lm  Tmin  values  of  Eq.  (41). 

Obtain  the  data  points  within  the  desired  BW  and  in  conjunction  with  the  data  from  Fig.  5(b) 
plots  of  Fig.  28  are  obtained. 


Step  4: 


Figure  28:  Frequency  Bounds  on  the  Prefilter  F(s). 

Utilizing  Fig.  28,  the  straight  line  Bode  technique  is  applied  and  with  the  condition 

limF(s)  =  1 

s-»0 


(42) 


for  a  step  forcing  function  satisfied  an  F(s)  is  synthesized  as  shown,  in  dashes,  that  lies  within 
the  upper  and  lower  plots  in  Fig.  28. 

Design  Procedure  for  the  QFT  Technique  for  M.P.  Plants 

Step  1 :  Synthesize  the  tracking  model  control  ratios  TR(s):  TR)J(s)  TRl(s) 

Step  2:  Synthesize  the  disturbance-rejection  model  control  ratios  TD(s) 

Step  3 :  Specify  J  FTI  plant  models  (define  the  plant  uncertainty) 

Step  4:  Obtain  templates  of  P(jec>i)  for  specified  <x>i  (describes  uncertainty)  within  the  BW 

Step  5 :  Select  a  nominal  plant  P0(s) 

Step  6:  Determine  the  stability  contour  (U-contour)  on  N.C. 

Step  7:  Determine  the  disturbance  bound  BD(jec>i)  on  N.C. 

Step  8:  Determine  the  tracking  bound  BR(jo)i)  on  N.C. 

Step  9:  Determine  the  optimal  bounds  Fm  B0(ja>i)  vs.  for  values  of  CO; 

Step  10:  Synthesize  the  nominal  L0(jcc>i)  =  G(jeo)P0(jeo)  that  satisfies  all  the  bounds  and  the  stability 

contours  and  then  obtain  G(s)  =  F0/P0(s) 

Step  1 1 :  Synthesize  the  prefilter  F(s) 

Step  12:  Simulate  the  FTI  system  (J  time  responses) 

Step  13:  Simulate  with  the  nonlinearities 
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1.3  MISO  Discrete  Control  System 


A  MISO  QFT  digital  control  system  is  shown  in  Fig.  29.  The  z-domain  equations  that  describe  this  system 
are: 


Pz(z)  =  GzoP(z)  =  (l— z-QZ 


m 

s 


=  (l-z-i)Pe(z) 


(43) 


L(z)  =  GZ0P(z)Gi(z) 


Pe (s)  — 


P(s) 

s 


(44) 


Pe(z)  -  2 


P(s) 

s 


2  [Pe (s)] 


Y*(s)  T 

Figure  29:  A  MISO  Sampled-data  Control  System. 


For  a  unit  step  disturbance,  D(s)  =  1/s: 


Pe(s)  =  P(s)D(s)  =  P(s)/s 


(45) 


PD(z)  =  Z[P(s)D(s)]=Pe(z) 


Tr(z)  = 


F(z)L(z) 
1  +  L(z) 


=  F(z)T'(z)  L(z) 


and 


Td(z)  = 


PD(z) 

1  +  L(z) 


(46) 


where 


t;cz)  = 


L(z) 

1  +  L(z) 


L(z)F(z)  PD(z) 

Y(z)  =  ———  R(z)  + - —  =  Yr  (z)  +  Yd  (z)  =  Tr  (z)R(z)  +  YD  (z) 

1  +  L(z)  1  +  L(z) 


(47) 


1.4  Pseudo-Continuous-Time  (PCT)  System  of  a  Digital  Control  System 

Introduction  —  The  PCT,  a  digitization  (DIG)  technique,  allows  the  QFT  design  of  the  z-domain 
controller,  D(z),  in  the  s-domain  [Houpis,  et  al.  (1999)].  That  is,  the  synthesized  s-domain  compensator 
(controller)  D(s)  is  transformed  into  the  z-domain  by  use  of  the  Tustin  transformation,  a  bilinear 
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transformation,  to  obtain  D(z).  The  advantage  of  this  approach,  especially  when  the  plant  is  m.p., 
eliminates  dealing  with  a  n.m.p.  plant  when  doing  the  analysis  and  design  in  the  w -plane  and  the 
associated  problem  in  satisfying  the  w'-domain  stability  bounds.  Note:  in  Fig.  29  D(z)  =  Gi(z). 

Another  important  advantage  is,  with  sampling  times  T  <  0.01  s  that  are  know  available,  the  waiping 
effect  that  is  discussed  later  on,  is  minimal  thus  the  s-domain  characteristics  are  transformed  into  the 
z-domain  assuring  that  the  robust  QFT  design  is  achieved  in  the  z-domain.  The  transformation  into  the 
s-domain  enables  the  use  of  the  MISO  QFT  analog  design  technique  to  be  used,  with  minor  exceptions, 
to  design  the  controller  D(s).  If  the  s-domain  simulations  satisfy  the  desired  P.S.  then  by  use  of  the  bilinear 
transformation  the  controller  G(z)  is  obtained.  With  this  synthesized  z-domain  controller  a  discrete-time 
domain  simulation  is  perfonned  to  verify  the  goodness  of  the  QFT  design. 

Bilinear  Transformations  —  The  Tustin  transformation,  a  bilinear  transformation,  transforms  an 
s-domain  function  into  the  z-domain  or  vice-versa.  A  function  of  z  is  substituted  for  sq  in  implementing 
the  Tustin  transformation.  This  function  is: 


s 


q 


2  1  -  z 


VT  \  +  lTl  j 


(48) 


The  advantages  of  the  Tustin  algorithm  are:  it  is  easy  to  implement,  the  accuracy  of  the  Tustin  z-domain 
transfer  function  is  good  compared  with  the  exact  z-domain  transfer  function,  and  the  accuracy  increases 
as  the  sampling  frequency  fs  increases  or  the  sampling  time  T  decreases.  The  Tustin  bilinear 
transformation  for  q  =  1  is  defined  as: 


g  =  2  1-z-1  _  2  z-1 
T  1  +  z"1  T  z  +  1 

and  can  be  equated  to  the  trapezoidal  integration  (s’1)  method. 

Mapping  Characteristics  -  For  the  s-  to  z-plane  mapping  Eq.  (49)  is  rearranged  to  yield 

1  +  sT/2 

z  = - 

1  -  sT/2 


(49) 


(50) 


Utilizing  the  mathematical  expression 


1  +  a 
1  -  a 


—  gj2tan  fi 


and  letting  s  =  jrasp  in  Eq.  (50)  the  following  expression  for  z  is  obtained: 


1  +  jcbspT/2 
1  —  jcbspT/2 


_  e  j2tan  _i  ®spT/2 


(51) 


Where  for  the  exact  ^transform  the  expression  for  z  is  z  =  e  J“spT ,  where  cosp  is  an  equivalent  s-plane 
frequency.  This  expression  is  substituted  into  Eq.  (51)  to  obtain: 


.  j®spT 


exp 


j2tan 


_  j  PbpT 


(52) 
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Equating  the  exponents  yields: 


or 


©spT 

2 


tan-1 


©spT 

2 


tan 


©spT  ©sp^ 


When  ©SpT/2  <  17°,  or  ~  0.30  rad,  then 

©sp  ~  ©sp 


(53) 


(54) 


(55) 


Thus,  in  the  frequency  domain  the  Tustin  approximation  is  good  for  small  values  of  ©spT/2,  that  is, 
minimal  waiping  is  achieved  for  the  imaginary  component. 

The  warping  of  the  real  pole  component  is  determined  as  follows.  Substitute  z  =  e0sPT  and  s  =  6sp  into 
Eq.  (50)  to  yield: 


e°spT 


1  +  d spT/2 

1  -  d SpT/2 


(56) 


Replacing  e  °spt  by  its  exponential  series  and  dividing  the  numerator  of  the  right-hand  side  by  its 
denominator  in  Eq.  (56)  results  in  the  expression 


1  +  cSpT  + 


dspT 

1  -  d spT/2 


(57) 


If  |aspT  |  »  (aspT  f/2  (or  1  »  |aspT/2|)  and  1  »  |dspT/2  ,  then  minimal  warping  occurs  and 


1  A  1 

1  1  2 

r  SP 1 

K  |°sp  |  «  q- 

(58) 

When  Eqs.  (55)  and  (58)  are  satisfied  the  Tustin  approximation  in  the  s-domain  is  good  for  small 
magnitudes  of  real  and  imaginary  components  of  the  variable  s.  The  shaded  area  in  Fig.  30  represents  the 
allowable  location  of  the  dominant  poles  and  zeros  in  the  s  plane  for  a  good  Tustin  approximation.  That  is, 
only  the  dominant  poles  and  zeros  of  the  open-loop  transfer  function,  that  determine  the  achievement  of 
the  desired  P.S.,  need  be  in  the  shaded  area.  Due  to  the  properties  and  its  ease  of  use, 
the  Tustin  transformation  is  employed  for  the  digitization  (DIG)  technique.  Figure  31a  is  the  mapping  of 
cb  =  2(tan  eospT/2)/T  [see  Eq.  (54)]  and  Fig.  31b  illustrates  the  warping  effect  of  a  pole  (or  zero)  when  the 
approximations  are  not  satisfied. 
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Figure  30:  Allowable  Location  (shaded  area)  of  Dominant  Poles  and 
Zeros  in  S-plane  for  a  good  Tustin  Approximation. 


Figure  31:  Map  of  cb  =  2(tan  rospT/2)/T.  (a)  Plot  of  Eq.  (52)  and  (b)  Warping  Effect. 


Characteristic  of  a  Bilinear  Transformation  -  In  general  the  bilinear  transformation  transforms  an 
unequal-order  transfer  function  (ns  ^  ws)  in  the  s-domain  into  one  for  which  the  order  of  the  numerator  is 
equal  to  the  order  of  its  denominator  (nz  =  wz)  in  the  z-domain.  This  characteristic  must  be  kept  in  mind 
when  synthesizing  G(s)  and  F(s). 

Introduction  to  PCT  System:  A  DIG  Technique  —  A  satisfactory  pseudo-continuous-time  (PCT) 
model  is  obtained  for  the  sampled-data  (S-D)  system  shown  in  Fig.  32.  The  PCT  approximation  of  the 
sampler  and  the  ZOH  units  of  Fig.  32  are  shown  in  Fig.  33b.  The  multiplier  1/T  attenuates  the 
fundamental  frequency  of  the  sampled  signal  and  all  of  its  harmonics.  These  approximations  are 
represented  by  the  linear  continuous-time  unit  Ga(s),  as  shown  in  Fig.  33c.  The  dominant  poles  and  zeros 
of  the  PCT  model  should  lie  in  the  shaded  area  of  Fig.  30  for  a  high  level  of  correlation  with  the 
S-D  system.  The  frequency  component  E  (jco)  of  the  sampled  signal  e*(t)  represents  E(]eo)  of  the 
continuous-time  signal  e(t)  multiplied  by  1/T,  where  all  the  side-bands  of  e*(t),  are  also  multiplied  by  1/T, 
are  attenuated  due  to  the  low-pass  filtering  characteristics  of  a  S-D  system. 
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Figure  32:  The  Uncompensated  Sampled-data  Control  System. 


Sampler  ZOH 

approximation  approximation 


(a)  (b) 


R(s) 


(c) 


Figure  33:  (a)  Sampler  and  ZOH;  (b)  Approximations  of  the  Sampler  and  ZOH; 
(c)  the  Approximate  Continuous-time  Control  System  Equivalent  of  Fig.  32. 


Using  the  first  Pade'  approximation,  the  transfer  function  of  the  ZOH  device,  when  the  value  of  T  is  small 
enough,  is  approximated  as  follows: 


GZo(s)  =  - - —  ~  2T  =Gpa(s)  (59) 

s  Ts  +  2 

where  the  transfer  function  Gpa(s)  is  used  to  replace  Gz0(s)  as  shown  in  Fig.33b  and  c.  Thus,  the  sampler 
and  ZOH  units  of  the  S-D  system  are  approximated  in  the  PCT  system  of  Fig.  33c  by  the  transfer  function: 

GA(s)  =  ~Gpa(s)  =  —  —  —  (60) 

T  Ts  +  2 


Since  Eq.  (60)  satisfies  the  condition  limT^)0GA(s)  =  1  it  is  an  accurate  PCT  representation  of  the  sampler 
and  ZOH  units.  That  is,  it  satisfies  the  requirement  that  as  T  — >  0  the  output  of  Ga(s)  must  equal  its  input 
and  in  the  frequency  domain  as  eos  — >  oo  (T  — >  0)  the  primary  strip  becomes  the  entire  frequency-spectrum 
domain,  the  representation  for  the  continuous-time  system. 


A  Simple  PCT  Example  —  As  an  illustration  of  the  effect  of  the  value  of  T  on  the  validity  of  the  results 
consider  the  S-D  closed-loop  control  system  of  Fig.  32  where 


Gx(s) 


Kx 

s(s  +  l)(s  +  5) 


(61) 


The  closed-loop  system  performance  is  determined  for  three  values  of  T  in  both  the  s-  and  z-domains 
[DIG  technique  vs  the  direct  (DIR)  technique  (the  z-plane  analysis)].  For  each  value  of  T  and  by  use  of  the 
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root-locus  technique  the  dominant  closed-loop  poles  are  determined  for  a  Q=  0.45.  Table  3  presents  the 
required  value  of  Kx  and  the  time-response  characteristics  for  each  value  of  T.  Note,  (a)  for  T  <  0. 1  s 
a  high  level  of  correlation  exists  between  the  DIG  and  DIR  models,  and  (b)  for  T  <  1  s  still  a  relatively  good 
correlation  exists.  (The  designer  needs  to  specify,  for  a  given  application,  what  is  considered  to  be  a  “good 
correlation.”) 


Table  3:  Analysis  of  a  PCT  System  representing  a  Sampled-data  Control  System  for  £  =  0.45 


Method 

T,s 

Plane 

Mp 

TP,s 

Ts,s 

DIR 

0.01 

z 

1.202 

4.16 

9.53 

DIG 

s 

1.206 

4.11 

9.478 

DIR 

0.01 

z 

1.202 

4.2+ 

9.8+ 

DIG 

s 

1.203 

4.33' 

9.90+ 

DIR 

1.0 

z 

1.199 

6 

13-14 

DIG 

s 

1.200 

6.18 

13.76 

The  MISO  PCT  QFT  Approach  —  The  PCT  System  of  the  digital  control  system  of  Fig.  29  is  shown  in 
Fig.  34a.  Note  that  the  sampler  sampling  the  forcing  function  and  the  sampler  in  the  system ’s  output  y(t) 
are  replaced  by  a  factor  of  1/T.  The  diagram  in  Fig.  34a  is  simplified  to  the  one  shown  in  Fig.  34b  and  is 
the  structure  that  is  used  for  the  QFT  design.  The  QFT  design  technique  can  now  be  applied  to  this  PCT 
system. 


PCT  Design  Summary  -  Once  a  satisfactory  Dc(.v)  [=  Gi(s)/T]  (see  Fig.  34b),  whose  order  of  the 
numerator  and  denominator  are,  respectively,  ws  and  ns,  has  been  achieved,  then  (a)  if  ns  is  equal  to  or 
greater  than  ws  +  2  then  the  exact  ^transform  is  applied  to  obtain  the  discrete  controller  Dc(z); 
or  (b)  if  ws  <  ns  <  ws  +  2  then  ns  -  ws  non-dominant  zeros  need  to  be  added  to  Dc(s)  before  applying  the 
Tustin  transformation  to  obtain  [Dc(z)]Tu.  As  a  final  check,  before  simulating  the  discrete  design,  obtain  the 
Bode  plots  of  Dc(s)|s  =  J(0  and  Dc(z)  |z  =  ejffl  .  If  the  plots  essentially  lie  on  top  of  one  another,  within  the 

desired  BW,  then  the  discrete-time  system  response  characteristics  are  essentially  the  same  as  those  for  the 
PCT  system  response.  If  the  plots  differ  appreciably  this  implies  that  warping  has  occurred  and  the  desired 
discrete-time  response  characteristics  may  not  be  achieved  (depending  on  the  degree  of  the  warping). 
If  the  warping  is  not  negligible  a  smaller  value  of  T  needs  to  be  selected  if  allowable. 


d(t),  D(s) 


(a) 


D(s) 


(b) 

Figure  34:  The  PCT  Equivalents  of  Fig.  29. 
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Accuracy  of  the  PCT  Design  —  An  analysis  of  the  characteristic  equation  (^(z)  =  1  +  Lt  [see  Houpis, 
et  ah,  (1999)]  reveals  the  following:  depending  on  the  value  of  T  and  the  plant  parameter  uncertainty  the 
pole-zero  configuration  in  the  vicinity  of  the  -1  +  jO  point  in  the  z-plane,  see  Fig.  35,  for  1  or  more  of  the 
J  plants  can  result  in  an  unstable  response.  Analyzing  the  characteristic  equation  Qv(z)  for  all  J  LTI  plants 
and  by  an  analysis  of  the  root-locus,  see  Fig.  36,  can  reveal  an  unstable  case(s).  This  analysis  revels  that 
by  a  slight  relocation  of  1  or  more  G(z)  nondominant  poles  in  the  vicinity  of  the  -1  +  jO  point  toward  the 
origin  may  ensure  a  stable  system  for  all  J  plants  and  maintain  the  desired  loop  shaping  characteristics  of 
G(z). 


X  -  Plant  poles  A  -  Controller  poles 

O  -  Plant  zeros  •  -  Controller  zeros 


Figure  35:  A  Root-locus  Pole-zero  Location  resulting  in  an  Unstable  System  Pole. 


Im  z 

A 


z-Plane 


-  Pole  of  G,(z)— Pg  (z) 


-1  + 


© 


■  Unit  Circle 


i  Re  z  (&) 

-4 — > 


i  +jo 


v_ 


Zero  of  Pe.(z)  =  zp 


(b) 


(c) 


Figure  36:  Analysis  of  Root-locus  Plot  of  L(z)  =  Gi(z)Pe,(z)  =  -1:  (a)  Pole-zero  Locations 
determining  System  Stability;  (b)  Location  yielding  a  Stable  System,  and  (c)  an  Unstable  System. 
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Accuracy  —  The  CAD  package  that  is  used  determines  the  degree  of  accuracy  of  the  calculations  and 
simulations.  The  smaller  the  value  of  T,  the  greater  the  accuracy  that  can  be  achieved  and  it  can  be 
enhanced  by  simulating  G(z)  and  F(z)  as  cascaded  functions,  that  is: 


G(z)  =  Gi(z)G2(z)-Gg(Z) 

(62) 

F(z)  =  F|  (z)F2  (z)  •  ■  •  Ff  (z) 

(63) 

2.0  MIMO  QFT  FUNDAMENTALS 

2.1  Introduction 

Part  2  discusses  how  the  QFT  technique  of  Part  1,  as  applied  to  MISO  control  systems,  can  be  applied  to 
the  design  of  MIMO  control  systems. 

2.2  System  Description 

Figure  37  represents  an  mxm  MIMO  closed-loop  system  where  F,  G,  and  P  are  mxm  matrices. 
The  notation  3{P}  denotes  the  set  of  J  matrices  (plant  parameter  uncertainty).  There  are  m2  closed-loop 
system  transfer  functions  (transmissions)  ty.  The  system  transmission  matrix  T  =  {ty}  relates  the  outputs 
y;  to  the  inputs  q,  that  is,  y;  =  typ. 


Figure  37:  MIMO  Control  Structure  for  a  Three-by-Three  System  (3x3). 


2.3  System  Analysis 

As  discussed  in  Part  1 ,  due  to  the  low  sampling  times  now  available,  the  controllers  for  a  digital  control 
systems  can  be  designed  by  use  of  the  PCT  approach.  Thus,  Part  2  restricts  the  analysis  to  the  s-domain. 
In  a  quantitative  problem  statement  the  tolerances  bounds  on  each  ty  specify  a  set  of  m2  acceptable  regions 
Ty  thus  the  design  sets  are: 

tij  £=  Ty  &  ^  {Ty} 

2.4  Introduction  to  Compensation 

The  system  control  ratio  relating  r  to  y  (see  Fig.  38)  is: 

T  =  [I  +  PG]  ‘  PGF  (64) 

where  the  ty  expressions  from  Eq.  (64)  are  complex  and  are  not  suitable  for  analysis.  Thus,  the  QFT  design 
procedure  systematizes  and  simplifies  the  achievement  of  a  satisfactory  system  design  for  the  entire  range 
of  plant  parameter  uncertainty.  In  order  to  readily  apply  the  QFT  design  technique,  another  mathematical 
system  description  is  used. 
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rl 


rm 


yi 


ym 


Figure  38:  A  MIMO  Feedback  Structure. 


2.5  Derivation  of  mxm  MISO  System  Equivalents  of  a  MIMO  System 

The  G,  F,  P,  and  P"1  matrices  are: 


gl 

0  ■ 

.  0~ 

"fn 

fl2  ' 

■  flm 

0 

g2  • 

.  0 

fni 

f22  ' 

•  f2m 

F  = 

0 

0 

•  §m_ 

/ml 

("m2  . 

f 

•  imm_ 

Pll  Pi 2  -  Plm 
P21  P22  •••  P2m 


Pml  Pm2  Pmm 


*  *  * 

Pll  Pl2  -  Pim 

*  *  * 

P21  P22  P2m 


P 


-1 


LPml  Pm2 


(66) 


(67) 


Note  that  only  a  diagonal  G  is  considered  for  this  paper  although  a  nondiagonal  G  allows  more  design 
flexibility.  The  m2  effective  plant  transfer  functions  are  based  upon  defining: 


q^i/p^ 


det[P] 
AdjyP 


(68) 
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where  the  requirement  det  P  be  m.p.  in  order  for  q,,  be  m.p.  The  Q  matrix  is: 


‘-111  *ll2  •”Cllm 

q2i  q22  ••■q2m 


Clm I  9  m2  Cln 


17P1 1 

1 7  P12 

-1/p* 

17  P21 

1/ p22 

—  IV  p2 

2m 


17  Pml  17  Pm2  -17Pr 


(69) 


The  inverse  plant  matrix,  P"1,  is  partitioned  to  the  form: 

P'1  =[Py]=[i/qij]  =  a+b  (70) 

where  A  is  the  diagonal  part  and  B  is  the  balance  of  P  1 . 

Thus 

hi=  1/qii  =  P*; ,  bji=  0,  and  for  i  *  j  by  =  1  /qy  =  pt 
Premultipling  Eq.  (64)  by  [I  +  PG]  yields: 

[I  +  PG]T  =  PGF 

Multiplying  both  sides  of  this  equation  by  P  *  yields: 

[P'1  +  G]T  =  GF  (71) 

where  P  is  nonsingular.  Using  Eq.  (70)  (G  diagonal)  Eq.  (71)  is  rearranged  to  the  form: 

T  =  [A  +  G]  '[GF  -  BT]  (72) 

Equation  (72)  defines  the  desired  Schauder  fixed  point  mapping  and  where  each  of  the  m2  matrix  elements 
on  the  right  side  of  this  equation  represents  a  M1SO  problem.  Thus,  the  design  of  each  M1SO  system  will 
yield  a  satisfactory  MIMO  design. 

Schauder’s  fixed  point  mapping  theorem  is  described  by  defining  a  mapping,  for  a  unit  impulse  forcing 
function,  as  follows: 

Y(T0  =  [A  +  G]“'  [GF  -  BT,]  =  T,  (73) 

where  each  T;  and  T,  are  from  the  acceptable  set  J.  If  the  mapping,  see  Fig.  39,  has  a  fixed  point, 
i.e.,  Tj,Tj  e  T  such  that  Y(Ty)  =  Tj,  then  a  robust  solution  has  been  achieved.  For  a  3x3  case,  for  a  unit 
impulse  input,  Eq.  (73)  yields  the  output: 


r  \ 

fill 

gifn- 

t21  .  t3l 

Ml  1  _L 

' 

1+§iqn 

V  9l2  fil3  y 
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Enclosure  of  all  acceptable  Te3 


Figure  39:  Schauder  Fixed  Point  Mapping. 


Based  upon  the  derivation  of  all  the  yy  expressions  from  Eq.  (73)  yields  the  effective  M1SO  loops 
(boxed  area  of  Fig.  40)  for  a  2x2  system  and  the  nine  effective  loops  resulting  for  a  3x3  system. 
The  control  ratios  for  the  desired  tracking  inputs  q  by  the  corresponding  outputs  y ,,  for  each  feedback  loop 
of  Eq.  (73),  for  unit  impulse  inputs,  have  the  form 

Yij =  wn  (vy  +  cy)  (75) 

where 


Wj,  =  q„/(l  +  g;qi)  and  Vy  =  g.f. 


-i  -i  -i 

Figure  40:  Effective  MISO  Loops  Two-by-Two  (boxed  in  loops) 
and  Three-by-Three  (all  nine  loops). 


The  interaction  between  the  loops  has  the  form 


c..  =-Z 
y  k*i 


Mk 


k  =  1,2,..., m 


(76) 


that  appears  as  a  “ cross-coupling ”  (disturbance)  input  in  each  of  the  feedback  loops. 


Equation  (75)  represents  the  control  ratio,  for  a  unit  impulse  forcing  function,  of  the  ith  MISO  loop  where 
the  transfer  function  wavy  relates  the  ith  output  to  ith  “desired”  input  q  and  the  transfer  function  waCy  relates 
the  ith  output  to  the  jth  “cross-coupling  effect”  input  cy. 
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Transfer  Functions  —  The  term  w^vy  relates  the  “desired”  ith  output  to  the  jth  input  r,  (see  Fig.  40) 
and  where  the  term  wuCy  relates  the  i“  output  to  the  jth  “cross-coupling”  input  cy.  Thus,  the  outputs 
[see  Eq.  (75)]  are  now  expressed  as 


Yij 

=  (yij)rj  +  (yij)cy  -  yrj  +  y,;j 

(77) 

or,  based  on  a  unit  impulse  input: 

Tracking 

> 

$ 

ii 

ii 

t 

(78) 

Cross-coupling 

lcii-yci,_wiiCij 

(79) 

Object  of  Design  —  The  design  objective  is  that  each  loop  tracks  its  desired  input  while  minimizing  the 
outputs  due  to  the  cross-coupling  inputs.  Note  that  each  MISO  system’s  cross-coupling  effect  input  is  a 
function  of  all  other  loop  outputs. 

MISO  Control  Ratios  —  To  achieve  the  desired  MISO  control  ratios  for  the  nine  MISO  structures  of 
Fig.  40  requires  that  the  ty  must  be  a  member  of  the  acceptable  ty  e  zy  and  the  g;  and  fy  must  be  chosen  to 
ensure  the  condition  on  ty  is  satisfied.  To  satisfy  these  constraints  constitutes  nine  MISO  design  problems. 
If  all  MISO  problems  are  solved  a  unique  solution  is  achieved  for  Fig.  37. 

2.6  Performance  Bounds 

The  developments  in  this  paper  are  based  on  the  utilization  of  diagonal  G  &  F  matrices,  i.e.: 

fij  =  gy  =  0  for  i*j 

In  Fig.  40  the  ty  terms,  for  i  ^  j,  involve  only  the  cross-coupling  (disturbance)  responses: 

yy  =  wyCij 

The  ty  terms,  for  i  =  j,  [see  Eqs.  (77)-(79)]  are  composed  of  a  desired  tracking  term  tr  and  an  unwanted 
cross-coupling  term  tc.  The  desired  tracking  specifications  are  given  by  the  upper  and  lower  bounds  of 
Fig.  41.  Whereas  the  cross-coupling  specification  for  all  MISO  loops  are  given  by  only  an  upper  bound. 


Figure  41:  Upper  and  Lower  Tracking  Bounds. 


Bounds  —  The  required  specified  P.S.  bounds  are  shown  in  Fig.  41  and  42.  The  upper  bounds  bii, 
within  the  BW,  is 

rCii  =  bii  -  bii  (80) 
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where  rqi  represents  the  maximum  portion  allocated  for  cross-coupling  rejection,  where  by  represents 
the  upper  bound  for  the  tracking  portion  of  t,;,  and  where  a',  represents  the  lower  bound  for  the  tracking 
portion  of  f,. 


Figure  42:  Allocation  for  Tracking  and  Cross-coupling  Specifications  for  the  tn  Responses. 

Modified  Tracking  Specifications  —  Based  upon  Eqs.  (78)  -  (80)  the  specifications  for  the  |tli(ja))| 
responses  of  Fig.  41  are  modified  as  shown  in  Fig.  42.  The  modification  is  based  upon  the  allocation  of 
SrOcoO  [see  Fig.  5b]  and  the  allocated  Tc;i  for  cross-coupling  effect  specification.  Based  upon  the  specified 
BW  above  which  the  output  sensitivity  is  ignored,  synthesize  the  G  and  F  such  that  for  all  P  s  IP: 

a;l=u,11)Lsitrlli^(!-rll)u-b:.  (8i) 

A  finite  oil,  is  recommended  since  in  strictly  proper  systems  feedback  is  not  effective  in  the  high  frequency 
range. 


2.7  Cross-Coupling  (Disturbance)  Specification  for  i  *  j 

Only  an  upper  bound  is  required  for  the  disturbance  specification,  i.e.: 

Itcijl  *  Ibijl  (82) 

Thus,  the  synthesis  of  G(s)  and  F(s)  must  satisfy  both  Eqs.  (81)  and  (82). 


2.8  Example:  A  2x2  MIMO  System 

The  ty  control  ratios  for  this  2x2  MIMO  system  example,  from  Eq.  (75)  for  unit  impulse  inputs,  are 

„  f  t2i  „  f  til 

gjlll" §2  121" 

c  4l2  421 

For  input  ri :  tn= — ; -  t2i  = - 


(83a) 


gifn-^k 
4 1 2 

1 


§2f22- 


tl3 

421 


1 

<122 


+  §2 


For  input  r2  :  ti2 


t22 


(83b) 


Multiplying  the  tn  and  t,2  equations  by  qn  and  the  t2i  and  t22  equations  by  q22  in  Eq.  (83),  respectively, 
yield  the  equations  shown  in  Fig.  43.  Associated  with  each  equation,  in  this  figure,  are  their  corresponding 


SFG. 


Figure  43:  A  2x2  MISO  Structures  and  their  Corresponding  t,j  Equations. 


Tracking  Bounds  —  The  ii  tracking  bounds  are  determined  in  the  same  manner  as  for  a  MISO  system; 
i.e.,  by  the  use  of  templates  for  the  ii  loop  plant  the  value  of  Axrii,  shown  in  Fig.  42,  are  used  to  satisfy  the 
constraints  ofEq.  (81). 

Cross-coupling  Bounds  —  From  Eqs.  (84)  and  (86),  in  Fig.  43,  the  following  cross-coupling  transfer 
functions  are  obtained: 


Cl  1 


cll% 


1  +  L 


-rcij“rcn 


(88) 


-12 


C12% 


1  +  L 


-bjj  =  t*!2 


(89) 


Substituting  for  cn  and  Ci2  into  Eqs.  (88)  and  (89),  then  replacing  t21  and  t22  by  their  respective  upper 
bound  values  b2i  and  b22,  and  then  rearranging  these  equations  will  yield  the  following  equations: 
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fil2 
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fil’ 
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(90) 


(91) 


In  order  to  use  the  N.C.  it  is  necessary  to  substitute  into  these  equations  Li  =  Ml \  which  results  in  the 
following  equations: 
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\  +  ix 


<MU 


\  +  £x 


<M12 


(92) 


These  equations  permit  the  determination  of  the  Mm  value  for  each  value  of  co;  over  the  set  J.  For  example: 
since  Li  =  1/Q  the  reciprocal  of  these  Mm(ja>j)  values,  based  on  lx,  yield  the  values,  based  on  L[  of  the 
corresponding  M-contours  or  the  cross-coupling  bounds  Bcl(jcOi)  for  co  =  w,  on  the  N.C.  In  order  to  use  a 
CAD  package,  Eqs.  (90)  and  (91)  are  rearranged  to 


|l  +  Li|> 


Mu 

9l2 


b21 

bll 


max  over  all  J  plants 


|l  +  Li 


> 


9n 

b22 

9l2 

b12 

max  over  all  J  plants 


Thus,  the  largest  magnitude  of  |qn  /qi2|  over  the  set  J  is  readily  determined  by  use  of  the  CAD  package. 
This  magnitude  is  inserted  into  Eqs.  (90)  and  (91)  to  yield  the  corresponding  cross-coupling  bounds  on  the 


N.C. 


Optimal  or  Composite  Bounds  —  The  points  on  the  composite  bound,  for  a  given  value  of  co;,  and  for  a 
given  row  of  MISO  loops  of  Fig.  40  is  the  composite  of  one  tracking  bound  and  m  cross-coupling  bounds 
for  this  value  of  C0i.  That  is,  select  the  largest  dB  value,  for  a  given  phase  angle  on  the  N.C.,  from  all  the 
tracking  and  cross-coupling  bounds  for  these  loops  at  this  value  of  frequency  to  be  the  point  on  the 
optimal  or  composite  bound.  The  M1MO  QFT  CAD  package  developed  by  Sating,  (1992)  is  designed  to 
perform  this  determination. 


2.9  Weighting  Matrix 

For  control  systems  whose  plant  matrix  Pb  has  more  inputs  V,,  than  outputs  Ym  a  weighting  matrix 
W  =  {Wij},  see  Fig.  44,  is  utilized  to  obtain  an  effective  square  mxm  plant  matrix  Pe,  i.e., 

Pe  =  PbW 


where  Pb  is  an  mxf  matrix  and  W  is  an  txm  weighting  matrix. 

r - 

u, - 


um 
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v, 
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p. 


Y, 


>  Y„, 


Figure  44:  An  mxm  Effective  Plant  Pe(s). 


2.10  QFT  Design  Methods 

There  are  two  QFT  design  methods  Method  1  and  Method  2.  Method  1  involves  synthesizing  the  loop 
transmission  L,  and  prefilter  fu  which  are  independent  of  the  previous  synthesized  loop  transmission  and 
prefilter  functions.  Whereas,  Method  2  involves  substituting  the  synthesized  g;  and  fb  of  the  first  (or  prior) 
MISO  loop(s)  that  is  (are)  designed  into  the  equations  that  described  the  remaining  loops  to  be  designed. 
At  the  onset  a  decision  needs  to  be  made  as  to  the  order  that  the  L,  functions  are  to  be  synthesized. 
Generally  loop  i  is  chosen  on  the  basis  of  the  performance  requirements.  For  example:  the  loop  i  having 
the  smallest  value  of  c»(|)  is  chosen  as  the  first  loop  to  be  designed  then  followed  by  designing  the  loop 
having  the  next  smallest  of  as  second  loop,  etc.  This  is  an  important  requirement  for  Method  2. 
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Method  1  —  Method  1  involves  over  design  (worst  case  scenario).  In  determining  the  My.  values  of 
Eqs.  (90)  and  (91)  (2  x  2  case),  for  each  value  of  co,,  obtain  the  smallest  magnitude  of  Iqn/qnl  (or  the  largest 
value  of  |qn/qi2|)  in  Eqs.  (90)  and  (91)  over  the  entire  J  LT1  plant  set.  These  values  are  utilized  to 
determine  the  bounds.  This  method  requires  that  diagonal  dominance  condition  be  met.  If  this  condition  is 
not  satisfied  Method  2  needs  to  be  used. 


Method  2  —  This  method  involves  selecting  the  loop  i  having  the  smallest  phase  margin  frequency 
requirement  to  be  designed  first.  For  example,  by  selecting  loop  i  =  1  results  in  the  synthesis  of  Gi  and  fn- 
These  are  now  known  LTI  functions  which  are  used  to  define  the  loop  i  =  2  effective  control  ratio  transfer 
functions.  That  is,  substituting  Eq.  (84)  [tn]  into  Eq.  (85)  [t2i]  and  rearranging  results  in  a  new  expression 
for  t2i  in  terms  of  gi  and  fn,  i.e.: 


t 


21 


fnLiq22e 
q2i(l  +  Lt) 
1  +  g2922e 


(93) 


The  effective  loop  2  transfer  function  is,  by  definition: 

q220  +  Li)  qnT22 

q„  = - - -  where  y12  = - 

1  +  (Li  -  Yi2)  qi2q2i 


Repeating  a  similar  procedure  the  expression  for  t22  is: 

f22g2q22e 

t22  = 

1  +  g2q22e 


Equations  (93)-(95)  involve  known  fn  and  gi  which  reduces  the  over  design  of  loop  2. 


(94) 


(95) 


2.11  Synthesizing  the  Loop  Transmission  and  Prefilter  Functions 

Once  the  optimal  bound,  for  each  o>i  have  been  determined,  for  each  Lj  loop,  then  the  synthesis  procedures 
of  Section  1 .2  for  determining  the  loop  transmission  and  prefilters  transfer  functions  are  used  to  determine 
the  gi  and  the  f„  for  the  MIMO  system. 


3.0  SYSTEM  WITH  EXTERNAL  DISTURBANCE  INPUTS 

3.1  Introduction 

The  previous  sections  have  dealt  with  just  the  tracking  problem.  This  section  deals  with  the  regulator 
problem,  i.e.: 

Input  tracking  commands:  r(t)  =  0 

External  disturbance  input:  d(t)  =  (dk(t)}  k  =  1,2,. ..pc 

The  development  of  this  section  permits  the  design  of  a  robust  multivariable  control  system  having  both 
types  of  inputs.  The  air-to-air  refueling  problem  is  used  as  the  basis  for  the  development  of  the  regulator 
problem  [Trosen,  (1993].  An  automatic  control  of  the  receiving  aircraft  during  aerial  refueling  operations 
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is  most  beneficial.  The  results  are  applicable  to  any  regulator  and  tracking-regulator  problem 
[Houpis  et  al.,  (1999)]. 

3.2  Problem  Statement 

The  design  problem  is  to  maintain  a  precise  position  of  the  receiving  aircraft  (receiver)  relative  to  the 
tanker  in  presence  of  such  disturbances  as  wind  gusts,  and  in  changes  of  mass  and  moments  of  inertia 
(see  Fig.  45)  and  to  prevent  a  disconnect  of  the  refueling  boom  due  to  excessive  changes  in  position. 
The  six-degrees-of-freedom  (uncontrolled)  A/C  models  are  used  to  develop  the  required  state  space 
models.  A  set  of  J  =  16  bare  A/C  plants  are  used  to  account  for  the  uncertainty  of  the  C-135B  aircraft 
during  an  air-to-air  refueling.  See  Floupis  et  al.,  (1999)  for  more  details  used  in  the  development  of  these 
models. 


Hinge  Point  (Tanker) 


Figure  45:  Control  Problem  Geometry. 


3.3  Disturbance  Modeling 

Disturbance  models  are  generated  by  developing  the  augmented  state-space  models  of  the  aircraft  in  the 
presence  of  wind  gusts  and  fuel  transfer  inputs.  The  external  disturbance  inputs  are  represented  by  the 
vector  d  and  contain  three  disturbance  components.  These  components  are: 

—  Pitch  plane  wind  induced  disturbance  rpitCh 

—  Lateral  channel  wind  induced  disturbance  Tiat 

—  Refueling  disturbance  Trf 

The  total  disturbance  modeled  is  given  by: 

Td  =  Tpitchdpitch  +  Tlatdlat  +  Trfdrf  (96) 

The  state  and  output  equations  are: 

x  =  Ax  +  Bu  +  rd  and  y  =  Cx  (97) 


3.4  Plant  and  Disturbance  Matrices 

Based  on  zero  initial  conditions  from  Eq.  (97)  the  following  equations  are  obtained: 

sX  =  Ax-Bu-Td  (98) 
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x  =  [sI-A]  Bu  +  [sI-A]  Td 


(99) 


y  =  Cx  =  C[sI-A]  Bu  +  C[sl-A]  Td  =  P(s)u  +  Pd(s)d 


(100) 


where 

P(s)  =  C[sl  -  A]  'b  (101a) 

Pd(s)  =  C[sl  -  A]  'r  =  {Pdy}  (101b) 

and  where  the  plant  model  PF  is  partitioned  into  the  two  matrices  (see  Fig  46)  P(s)  and  Pd(s)  where 
P(s)  =  Pe(s)  for  a  square  plant  matrix  P(s)  and  the  matrix  Pd(s)  models  the  transmission  from  the  external 
disturbance  inputs  to  the  output  of  PF.  If  P(s)  is  not  a  square  matrix  then  a  weighting  matrix  W(s)  is  used 
to  yield  Pe(s)  =  P(s)W.  Equation  (100)  is  represented  in  Fig.  46. 


Figure  46:  QFT  Compensator  with  Output  External  Disturbance. 


3.5  MIMO  QFT  with  External  Output  Disturbance 

Figure  46  represents  an  mxm  closed-loop  system  in  which  F(s),  G(s),  P(s),  and  Pd(s)  are  mxm  matrices. 
Note  that  tP(s)  =  {P(s)|  and  (Pd(s)  =  {Pd(s)|  are  sets  of  matrices  due  to  the  plant  and  disturbance 
uncertainties,  respectively.  The  objective  is  to  find  a  suitable  mapping  that  permits  the  analysis  and 
synthesis  of  a  control  system  by  a  set  of  equivalent  M1SO  control  systems. 

Mathematical  Development  —  From  Fig.  46  the  following  equations  are  obtained: 

y(s)  =  P(s)u(s)  +  Pd  (s)  dext  (s)  u(s)  =  G(s)v(s)  v(s)  =  r(s)  -  y(s)  (102) 


For  the  regulator  case,  with  zero  tracking  input: 

r(t)  =  [o,0,of 


(103) 


From  Eqs.  (102)  and  (103),  where  henceforth  the  (s)  is  dropped  in  the  continuing  development,  obtain  the 
following  equations: 


which  yields 

Equation  (105)  is  rearranged  to  yield: 


v  =  -  y  u  =  -  Gy 

y  =  -  PGy  +  Pddext 

y  =  [l  +  PG}1Pdd 


ext 


(104) 


(105) 


(106) 
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Based  upon  unit  impulse  disturbance  inputs  for  dext  the  system  control  ratio  relating  dext  to  y  is 

Td  =  [l  +  PG]_1Pd 

Pre-multip lying  Eq.  (107)  by  [I  +  PG]  yields 

[l  +  PG]Td=  Pd 

Pre-multiplying  both  sides  of  Eq.  (108)  by  P  results  in 

[P-1  +  G]Td  =  P  % 


(107) 


(108) 


(109) 


Let 


Pi 


Pll 

Pl2  • 

•'  Plm 

P21 

P22  • 

•  •  P2m 

P  ml 

P^  • 

• '  Pmm 

The  m2  effective  plant  transfer  functions  are  formed  as 

1  det  P 

Oii  * 

J  Pij  adjPy 


(110) 


(111) 


The  Q  matrix  is  then  formed  as 


Q  = 

Oil  O12  ' 

921  922  ' 

'  9lm 

'  92m 

— 

'l/Pn  1/P n  • 
1/P21  !/P22  ' 

•  i/pL  " 

•  1/pL 

_0ml  9  m2  ' 

Omm  _ 

j/pim  i/p;2  • 

1  /P  mm  _ 

(112) 


where  P  =  [py],  P'1  =  [p*  ]  =  [1/qy],  and  Q  =  [qy]  =  [1/p*  ].  The  P  *  matrix  is  partition  as  follows: 


P'  =  [pJ]  =  [l/qy]  =  A  +  B  (113) 

where  A  is  the  diagonal  part  of  P  1  and  B  is  the  balance  of  P  1  Thus, 

A.ii  =  1/qii  =p*j ,  bn  =  0,  and  by  =  1/qy  =  p*.  for  i  *  j 

Substituting  Eq.  (113)  into  Eq.  (109)  with  G  diagonal,  results  in 

[A  +  B  +  G]Td  =  [A  +  B]Pd  (114) 


Rearranging  Eq.  (114)  produces 

Td  =  [  A  +  G  ?[  APd  +  BPd  -  BTd  ]  =  {  td|j}  (115) 


1  -40 


RTO-EN-SCI-142 


Quantitative  Feedback  Technique  (QFT):  Bridging  the  Gap 


This  equation  defines  the  desired  fixed  point  mapping,  where  each  of  the  m2  matrix  elements  on  the  right 
side  of  Eq.  (115)  are  interpreted  as  M1SO  problems.  The  theorem  defines  a  mapping  Y(Tdi) 

Y(Tdi)  =  [A  +  G]-J[APd  +  BPd  -  BTd.  ]  =  Tdj  (116) 

where  each  Td;  and  Tdj  is  from  the  acceptable  set  3d.  If  this  mapping  has  a  fixed  point,  i.e.,  Tdj),Tdj  s  3d 
such  that  Y(Tdi)  =  Tdj  then  a  robust  solution  has  been  achieved. 


Figure  47  shows  the  effective  MISO  loops  resulting  for  a  3x3  system.  Since  A  and  G  in  Eq.  (115) 
are  diagonal,  the  (1,1)  element  on  the  right  side  of  Eq.  (116)  for  the  3x3  case,  for  a  unit  impulse  input, 
provides  the  output 


(117) 


Figure  47:  Effective  MISO  Loops  Two-by-Two  (boxed  in  loops) 
and  Three-by-Three  (all  nine  loops). 


The  three  ratios  involving  the  elements  pd|  are  the  additional  terms  due  to  the  external  disturbances. 
The  remaining  terms  within  the  parentheses  are  the  same  terms  as  for  tracking.  Equation  (117)  corresponds 
precisely  to  the  first  structure  in  Fig.  47.  Similarly,  each  of  the  nine  structures  in  this  figure  corresponds  to 
one  of  the  elements  of  Y(Tdi)  of  Eq.  (116). 

The  control  ratios  for  the  external  disturbance  inputs  dext,  and  the  corresponding  outputs  y,  for  each 
feedback  loop  of  Eq.  (115)  have  the  form 

yii  =  wii(deij)  (ns) 


where  wd  =  qij/(l  +  g,q,)  and 


(dext  )ij  "  Cij  X 


Pd 


kj 


t. 


kj 


eii 


k=i  qik 


-  X 

k^'  4 ik 


(119) 


where  x  =  the  number  of  disturbance  inputs  and  m  =  the  dimension  of  the  square  M1MO  system. 
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Thus,  the  interaction  term,  Eq.  (119),  not  only  contains  the  cross-coupling  interaction  but  also  the  external 
disturbances,  i.e.: 


(dext  )ij 


k=r 


4ik 


and 


Cij  - 


k^ifiik 


(120) 


where  (dext)y  represents  the  external  disturbance  effects  and  c„  represents  the  cross-coupling  effects. 
The  additional  equations,  quantifying  both  the  external  disturbance  (dext)ij  and  the  internal  cross-coupling 
effects  cy,  are  derived  to  utilize  the  improved  method  QFT  design  technique.  These  equations  are  used  to 
define  the  disturbance  bounds  for  subsequent  loops  based  on  the  completed  design  of  a  single  loop. 


3.6  Tracker/Regulator  MIMO  Control  System 

A  MIMO  system,  that  involves  both  tracking  and  external  disturbance,  can  now  be  designed  by  either 
QFT  Method  1  or  Method  2.  In  Fig.  48  the  cy  terms  are  replaced  by  the  de]|  terms;  that  is,  Fig.  47  is 
modified  by  including  the  tracking  inputs  r,  ^  0.  The  modified  figure  now  represents  the  tracking/regulator 
MIMO  control  system.  The  optimal  bounds  are  now  a  combination  of  not  only  the  most  stringent  portion 
of  each  of  the  tracking  and  cross-coupling  effect  bounds  but  also  that  of  the  bounds  due  to  the  external 
disturbance  (dext)ij  of  Eq.  (120). 


Figure  48:  Effective  MISO  Loops  Two-by-Two  (boxed  in  loops)  and  Three-by-Three  (all  nine). 


4.0  NOW  THE  “PRACTICING  ENGINEER  TAKES  OVER” 

4.1  Introduction 

Scientific  Method  —  The  scientific  method  uses  mathematical  methods  to  gain  insights  into,  generalize, 
and  to  expand  the  state-of-the-art.  This  method  often  requires  proof  of  theorems,  corollaries,  and  lemmas. 
At  this  stage,  in  general,  the  researcher  (a)  is  not  concerned  whether  his/her  efforts  result  in  the  solution  of 
“real  world  problems”  which  are  nonlinear  and  (b)  applies  linear  analysis  and  synthesis  techniques  in  the 
research  effort.  The  scientific  method  is  necessary  for  the  researcher  to  be  able  “to  see  the  trees  from  the 
forest,”  and  thus  to  be  able  to  achieve  positive  results. 

Engineering  Method  —  The  engineering  method,  where  applicable,  is  a  method  where  the  engineer  takes 
over  and  applies  these  new  results  to  real  world  problems.  The  engineer  is  at  the  “ interface ”  of  the  real 
world,  and  with  the  body  of  knowledge  and  theoretical  results  available  in  the  technical  literature. 

Application  of  the  Scientific  Method  to  OFT  —  Professor  1.  M.  Horowitz  applied  the  scientific  method 
in  his  development  of  QFT  approach  to  the  engineering  design  of  robust  control  systems  [Horowitz,  et  al.. 
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(1973)].  Professors  F.  Bailey,  0.  D.  I.  Nwokah,  et  al.,  also  used  this  approach  to  enhance  the  mathematical 
rigor  of  QFT  and  to  help  the  engineer  to  bridge  the  “ interface ”  gap  in  applying  the  QFT.  This  body  of 
knowledge  must  be  coupled  with  the  “ body  of  engineering  knowledge ”  pertaining  to  the  application,  when 
dealing  with  nonlinear  systems  and  real  world  problems.  This  requires  that  the  engineer  have  a  good 
understanding  of  the  physical  characteristics  of  the  plant  to  be  controlled. 

“Bridging  of  the  Gap”  —  The  “ Bridging  of  the  Gap ”  [Floupis,  (2002)]  is  best  illustrated  by  the  following 
anonymous  saying: 

“In  THEORY  (Scientist) 

There  is  no  difference  between  THEORY  and  PRACTICE. 

In  PRACTICE  (Engineer) 

There  is  a  difference  between  PRACTICE  and  THEORY.” 

Conclusion  —  The  control  system  design  task  is  a  multi-stage  process  which  entails  many  steps,  say  from 
A  to  Z.  Mathematics  is  useful  in  taking  the  engineer  through  some  of  these  steps,  say  from  P  to  S. 
The  engineer  should  make  the  required  modelling  assumptions,  hypotheses,  and  simplifications  needed  to 
proceed  from  A  to  0,  so  that  the  mathematical  problem  is  tractable  and  the  existing  theory  can  be  applied. 
Finally,  the  steps  T  to  Z  entail  extensive  simulations  where  the  validity  of  the  model  is  verified, 
the  implementation  issues  are  addressed,  and  the  design  is  validated. 

4.2  Transparency  of  QFT 

There  are  essentially  four  elements  that  provide  the  “transparency”  of  the  QFT  technique  and  which 
enhance  its  ability  to  design  and  to  solve  real  world  problems  [Houpis,  et  al.,  (1999)].  These  elements  are: 

Size  of  Template  —  The  width  and  height  of  a  template,  see  Fig.  49,  play  an  important  role  in  Bridging 
the  Gap  at  the  onset  of  the  design  process.  That  is,  the  engineer  at  the  onset  of  the  design  process  can 
readily  determine  whether  a  fixed  compensator  G  can  be  synthesized  that  yields  the  desired  system 
performance.  If  only  the  template  height  is  the  problem  then  the  engineer  needs  to  employ  straight  gain 
scheduling.  When  control  effector  failures  need  to  be  accommodated,  the  width  of  the  template  can  be 
excessive  and  a  successful  design  is  not  possible  [Cacciatore,  (1995)].  For  this  case,  the  designer  needs  to 
either  reduce  the  percentage  of  failure  to  be  accommodated  or  to  eliminate  the  effector  failure  case(s), 
in  order  to  reduce  the  width  of  the  template.  Chapter  10  by  Houpis,  et  al.,  (1999)  provides  an  example  of 
how  an  engineer  utilized  the  knowledge  of  plant  parameters  and  the  performance  characteristics  and 
specifications,  to  select  a  set  of  J  plants  that  ensures  a  template  shape  which  will  represent  as  accurately  as 
possible  the  extent  of  plant  parameter  uncertainty. 


Figure  49:  A  Template  Representing  J  LTI  Plants  at  a  =  30  rad/sec. 
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Phase  Margin  Frequency  —  In  order  to  ensure  that  the  value  of  the  specified  phase  margin  frequency  co(|l 
is  not  exceeded  by  any  of  the  plants  in  the  set  then: 

(a)  When  all  J  plants  are  stable,  select  the  nominal  plant  p,  to  be  the  plant  lying  at  the  “top”  vertex  of  the 
template.  This  assumes  that  this  plant  always  lies,  for  all  template  frequencies,  at  the  top  of  the 
templates.  If  not,  next  best  use  the  template  for  co^,  to  select  the  “top”  plant  as  the  nominal. 

(b)  If  one  or  more  plants  in  the  set  9  are  unstable,  then  the  design  engineer  must  select  as  the  nominal 
plant  the  plant  P,  that  has  the  highest  degree  of  instability  (the  plant  whose  unstable  pole  lies  furthest 
to  the  right  in  the  s-plane)  as  the  “worst  case  plant.” 

This  rule  for  the  selection  of  the  nominal  plant  facilitates  the  achievement  of  the  specified  value  of  co^.  It  is 
advisable  that  one  of  the  templates,  for  each  loop,  be  obtained  at  the  respective  BW(L,)  frequency. 

Signal  Flow  Graph  —  The  signal  flow  graph  (SFG)  for  the  portion  of  Fig.  50  that  represents  Pe  is  helpful  in 
the  initial  selection  of  the  values  of  wy  of  the  weighting  matrix  W  [see  Eq.  (121)]  and  in  modifying  some 
of  these  values  during  the  simulation  phase  of  the  design  process. 

y  =  Peu  =  PWu  (121) 


R  Pe  Y 


Figure  50:  MIMO  QFT  Control  Structure  Block  Diagram. 

The  selection  of  values,  etc.,  can  be  further  enhanced  by  the  engineer’s  firm  understanding  of  the  interrelationship 
of  the  plant  outputs  with  the  inputs  to  the  W  matrix. 

Minimum  Order  Controller  (MOC)  G  —  The  implementation  of  an  m-order  controller  in  a  digital  flight 
control  system  (FCS)  results  in  m-time  delays.  Thus,  during  the  first  m-1  time  instants  the  input  has  a 
limited  effect  on  the  output,  for  the  output  is  partially  determined  by  the  m  initial  conditions.  Thus,  the 
control  action  is  delayed.  A  factor  to  be  considered  in  maintaining  low-order  controllers  is  that  on  board 
flight  control  computers  have  a  limited  capacity  due  to  other  non-control  related  computing  requirements. 
As  a  rule,  about  30%  of  the  computer  capacity  is  allocated  to  the  FCS  and  70%  for  non  FCS  requirements. 
To  achieve  this  desideratum  some  designers,  or  design  methods,  have  recourse  to  the  plant  P,  and  in-tum 
Pe  by  “doctoring”  or  “padding”  by  inserting  additional  poles  and/or  zeros  into  P.  This  results  in  an 
“augmented”  plant  matrix  Pes  which  these  designers  base  their  design  on  in  order  to  achieve  a  so-called 
“minimum- order”  G.  These  “additional  poles  and/or  zeros”  in  reality  are  poles  and  zeros  of  G. 

To  achieve  the  MOC  G,  during  the  QFT  loop  shaping  phase,  the  poles  and  zeros  of  the  nominal  plant 
qiio  are  used  in  synthesizing  a  satisfactory  loop  shaping  transfer  function  L;  .  Doing  so,  yields  the 
MOC  gi  =  Li  /qjj  Some  or  all  of  these  elements,  or  comparable  ones,  are  not  available  in  other, 
optimization  based,  multivariable  control  system  design  techniques.  This  minimizes  their  ability  to 
achieve,  in  a  relatively  “short  design  time,”  a  design  that  meets  all  the  performance  specifications  that  are 
specified  at  the  onset  of  the  initial  design  effort. 
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4.3  Body  of  Engineering  QFT  Knowledge 

Y ears  of  applying  the  QFT  robust  control  design  technique  to  real  world  nonlinear  problems  have  evolved 
the  following  Engineering  Rules  (E.R.).  These  rules  were  developed  by  Floupis,  et  al.  (1999),  and  the 
reader  is  referred  to  this  text  for  a  more  complete  presentation.  The  rules  are: 


E.R.1  Weighting  Matrix  —  The  weighting  matrix  W  =  {wy}  is  required  to  achieve  a  square  equivalent 
plant  matrix.  It  is  desired  to  know  at  the  onset,  if  possible,  to  achieve  m.p.  q^’s  by  the  proper  selection  of 
the  Wy  elements.  A  m.p.  qn  plant  is  most  desirable  for  it  allows  a  full  exploration  of  the  “benefits  of 
feedback,”  i.e.,  high  gain.  By  applying  the  Binet-Cauchy  theorem,  as  shown  by  Floupis,  et  al.,  (1999),  it  is 
possible  to  determine  if  m.p.  q^’s  are  possible.  It  may  be  desirable  to  obtain  complete  decoupling  for  the 
nominal  plant  case,  i.e., 


ediag 


Pll  0 

o  p22 

0  0 


PW 


(122) 


For  non-nominal  plants  complete  decoupling,  in  general,  will  not  occur  but  the  degree  of  decoupling  will 
have  been  enhanced.  Method  1  is  then  more  readily  applicable,  with  the  additional  benefit  of  reduced 
closed-loop  BW. 


E.R.2  n.m.p.  q^’s  —  For  manual  flight  control  systems,  if  a  RFIP  zero  is  “close”  to  the  origin  (see  Fig.  5 1) 
it  is  not  necessary  deleterious  since  the  pilot  will  input  a  new  command  before  the  effect  of  the  closed- 
loop  pole,  associated  with  this  zero,  is  noticeable.  Assume  that  the  unstable  pole  is  outside  the  closed  loop 
system’s  BW.  If  this  RHP  zero  is  “far  out”  to  the  right,  outside  the  BW  of  concern,  in  manual  control  it 
does  not  present  a  problem.  For  these  cases  a  satisfactory  QFT  design  may  be  achievable. 


Figure  51:  Right-half-Plane  Analysis. 

E.R.3  Templates  —  The  adage  “picture  is  worth  a  thousand  words”  applies  to  the  preliminary  task  of 
determining  if  a  robust  control  solution  exists.  If  theorems,  corollaries,  and/or  lemmas  obtained  by 
the  scientific  method,  reveal  that  no  loop  shaping  solution  exists  then  one  must  step  back  and  do  a 
“trade-off’  in  which  some  specifications  are  relaxed  in  order  to  achieve  a  solution,  or  one  must  be 
willing  to  live  with  a  degree  of  gain  scheduling.  A  graphical  analysis  can  reveal  the  following: 

(a)  The  maximum  template  height  is  too  large  forcing  one  or  more  bound  or  the  composite  bound  being 
violated.  The  engineer  needs  to  decide  if  gain  scheduling  is  required  and  if  it  is  feasible  in  order  to 
yield  a  design  that  satisfies  all  the  bounds. 
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(b)  The  analysis  may  reveal  that  the  templates  are  too  “wide”  (magnitude  of  the  phase  angle  width) 
thus  prohibiting  a  QFT  solution  or  a  solution  by  any  other  multivariable  design  technique.  This  is  true  for 
real-world  control  problems  that  involve  control  effector  failures.  In  these  design  problems,  generally, 
the  worst  failure  case  is  the  culprit  in  generating  large  “angle  width.”  Thus,  in  order  to  achieve  a 
solution  it  is  necessary  to  relax  the  requirement  that  the  “worst  failure  case”  be  accommodated.  It  is 
also  necessary  to  stipulate  what  failure  case  or  cases  a  successful  design  is  achievable.  In  determining 
the  “reasonable  failure  cases”  that  can  be  accommodated  by  robust  (not  adaptive)  control  one  must 
consider  if  10%,  ...  ,  80%  failure  still  permits  enough  control  authority!  This  percentage  of  failure  can 
only  be  determined  by  a  individual  who  is  knowledgeable  of  the  physical  plant  being  controlled. 
A  knowledge  of  the  plant  (application)  “is  king”  when  it  comes  to  the  design  of  a  feedback 
compensator  or  controller  for  the  given  plant. 

(c)  In  Section  10-4  of  Houpis,  et  al.  1999,  it  is  shown  how  the  engineer’s  insight  into  the  plant’s  physics 
helps  to  efficiently  determine  the  boundary  of  the  template  (where  the  parametric  uncertainty  is 
represented).  A  possible  method  of  reducing  the  size  of  the  templates  is  given  by  E.R.8. 

E.R.4  Design  Techniques  —  For  all  the  design  methods  the  P.S.  must  be  realistic  and  commensurate  with 
the  real  world  plant  being  controlled.  Situations  occurred  where  a  conclusion  was  reached  that  no 
acceptable  design  was  possible.  For  these  situations  when  one  “stepped  back”  and  asked  the  pertinent 
question  “was  something  demanded  that  this  plant  physically  cannot  deliver  regardless  of  the  control 
design  technique?”  It  was  determined  that  some  or  all  of  the  prescribed  specifications  were  unrealistic. 

E.R.5  QFT  Method  2  -  In  using  Method  2  it  is  known  at  the  onset  the  co^  of  the  succeeding  designed  loop 
is  larger  than  the  previously  designed  loops.  Thus,  arbitrarily  picking  the  wrong  order  of  the  loops  to  be 
designed  can  result  in  the  nonexistence  of  a  solution.  This  may  occur  if  the  solution  process  is  based  on 
satisfying  an  upper  limit  of  the  phase  margin  frequency  co(|),  for  each  loop.  The  proper  order  of  loops  to  be 
designed  by  Method  2  is  to  pick  the  loops  in  the  order  of  increasing  values  of  the  desired  oy,;  i.e.,  first  close 
loop  1,  then  loop  2,  etc  where  eo^  <  w{|n  <  co^3  <  ... . 

E.R.6  Minimum  order  G  —  To  ensure  the  smallest  possible  order  compensator/controller  it  is  necessary  to  start 
the  loop  shaping  process  by  using  the  loop’s  nominal  plant  as  the  starting  L0,  i.e.,  L0  =  L0i  =  qUo.  The  loop  shaping 
process  is  then  started  by  inserting  zeros  and  poles  successively  in  order  to  obtain  the  required  loop  shape, 
resulting  in: 


L0(s) 


L01(s)(s-zi)...(s-zw) 
(s  —  P ! )  •  •  •  (s  —  Pv) 


(123) 


Finally,  from  the  synthesized  loop  shaping  function  the  compensator  is  obtained,  i.e.,  gi  =  [Lo]i/qn0. 

E.R.7  Minimum  Compensator  Gain  —  To  minimize  the  effects  of  noise,  saturation,  etc.  need  to 
minimize  the  gain  required  in  each  loop  i  while  at  same  time  satisfy  the  P.S.  To  achieve  this  goal, 
a  designer,  with  a  good  understanding  of  the  N.C.  and  a  good  interactive  QFT  CAD,  can  use  his 
“engineering  talent”  to  make  use  of  the  “dips”  in  the  composite  Bol(jc0j)  -  See  Fig.  52.  Shaping  [L0]i 
to  pass  through  the  dips,  where  feasible,  ensures  achieving  the  minimum  compensator  gain  that  is 
realistically  possible. 
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Figure  52:  QFT  Stability  and  Composite  Bounds. 

E.R.8  Basic  mxm  Plant  P  Preconditioning  —  Using  the  unity  feedback  loops  for  the  mxm  M1MO  plant 
P  yields  an  mxm  preconditioned  plant  matrix  Pp.  The  templates  3Pp(jec>j),  in  general  can  be  smaller  in 
size  than  the  templates  3P(jcc>i).  Thus,  the  QFT  design  may  be  performed  utilizing  the  preconditioned 
matrix  Pp. 

E.R.9  Nominal  Plant  Determination  —  The  phase  margin  angle  y,  the  gain  margin,  and  o)(il  of  a  feedback 
control  system  can  be  readily  determined  using  the  N.C.  The  QFT  design  technique  affords  the  robust 
establishment  of  these  FOM.  By  choosing  the  nominal  plant  to  correspond  to  the  maximum  dB  plant  on  the 
co^i  =  BW(Lj)  template  ensures  that  co^,  <  BW(L;)  for  all  plants.  The  achievement  of  a  robustly  guaranteed 
gain  margin  is  accomplished  readily  by  selecting  a  nominal  plant  that  is  uniformly  the  maximum  dB  plant 
for  all  templates  at  the  top  of  the  template.  The  achievement  of  a  robustly  guaranteed  desired  y  is 
accomplished  readily  by  selecting  the  “left-most”  plant  on  all  the  templates  as  the  nominal  plant. 
This  ensures  that  the  desired  y  is  robustly  achieved. 

E.R.10  Simulation  Run  Time  —  The  goodness  of  the  design  of  the  “real-world”  manual  feedback 
control  problems  are  judged  on  a  pre-specified  planning  time  horizon  beyond  which  the  performance  is 
less  important  since  the  human  operator  will  inject  new  inputs  to  the  system.  For  example,  in  many 
manual  flight  control  scenarios  the  time  horizon  is  5  s  simulation  run  time. 

E.R.11  Asymptotic  Results  —  Asymptotic  results  by  mathematical  analysis  are  not  as  useful  as  they 
seem  to  be,  i.e.,  for  external  disturbance  yD(co)  =  0.  Consider  the  manual  control  disturbance  rejection 
case  where  fast  disturbance  attenuation  is  more  desirable  than  total  disturbance  rejection,  i.e.,  YD(oo)  =  0, 
which  entails  a  very  long  “settling  time.” 

E.R12  Controller  Implementation  —  To  minimize  the  computer  numerical  inaccuracy  due  to  the  software 
implementation  of  G(z),  it  is  necessary  to  obtain  the  software  implementation  by  writing  the  difference  equation 
for  each  block  in  Fig.  53. 
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I(s)- 


r  =  Mi) 

^c(S)  D(s) 


'O(s) 


I(w)- 


O(w) 


(a)  The  design  compensators 


I(z)- 


Gc(z)=^ 


O(z) 


I(z) 


Gj(z) 


(b)  z-domain  compensator 
Oj(z) 


»5 


S- 


Gj(z) 


O.(z)  =  O(z) 

- i - ^ 


(c)  z-domain  J  cascaded  controllers 

Figure  53:  s-  or  w-domain  to  z-domain  Bilinear  Transformation: 
Formulation  for  the  Implementation  of  the  G(z)  Controller. 


E.R.13  Non-ideal  Step  Function  for  Simulation  —  For  mathematical  analysis  of  LT1  systems  ideal  step 
functions  are  utilized.  In  the  real  world,  the  ideal  step  functions  do  not  exist.  Thus,  for  a  more  realistic  test 
in  determining  a  control  system’s  performance  through  simulation  a  ramped-up  step  function  is  used  as 
shown  in  Fig.  54.  Use  of  this  type  of  a  realistic  forcing  function  has  the  tendency  to  minimize  the 
saturating  aspects  of  the  system. 


Figure  54:  A  Ramped-up  Step  Function. 


4.4  Plant  Inversion 

Given  the  mxl  output  vector  y  and  an  mxl  input  vector  u  (see  Fig.  55)  the  LTI  plant  state  equation  is: 

D(s)y  =  N(s)u  (124) 

with  D(v)  =  [dij(.s)]  and  N(.v)  =  [n,j(.v)]  being  mxm  polynomial  matrices  in  s.  The  resulting  plant  matrix 
P  (y  =  Pu)  is 

(AdjD)N 

P  =  D-iN  =  - — - —  —  =  [p-  (s)]  (125) 

detD  J 
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Yi 

'  y  2 


Figure  55:  A  MIMO  Plant. 


The  inverse  plant  is  expressed  as: 

p  1  =  [Py]  =  [1/qij] 

It  can  be  obtained  from  either  Eq.  (124),  i.e.,  since  u  =  P  y, 

p-1  =  N-1D=(AdjN)D 
det  N 

or  from  the  state  equations,  i.e.: 

x  =  Ax  +  Bu  (a)  y  =  Cx  (b) 

that  describe  the  nth-order  plant,  where  the  matrices  A,  B,  C  are  nxn,  nxm,  and  mxn,  respectively,  and 
which,  in  turn,  yield  the  following  expressions: 

y  =  C[sI  -A]  'Bu 
P  =  C  [si  -  A]  -'B  =  jyj 

where  ny  and  d  are  polynomials  of  degree  l  and  n,  respectively,  in  s;  l  <  n  and  d  is  the  characteristic 
polynomial  of  A. 

For  the  numerical  calculation  of  P  *  it  is  much  better  to  use  Eq.  (127)  rather  than  Eq.  (130).  This  is  best 
illustrated  by  considering  a  2x2  (m  =  2)  plant  utilizing  the  second  approach.  Thus,  from  Eq.  (130): 

AdjW 

|dj  Adj  {ny}/d  d“Adj{ny} 

P  = - y  = - 2  = -  (131) 

det{n;j}/cf  det{ny}/cf  ddet{ny} 

In  general,  for  the  mxm  control  system: 


(129) 

(130) 


(126) 


(127) 


(128) 


dmAdj{niJ| 
d™  '  det  {ny} 


(132) 


Thus,  if  P  *  is  obtained  from  Eq.  (132),  rather  than  directly  from  Eq.  (127),  then  m-1  cancellations  of  the 
polynomials  from  the  numerator  in  Eq.  (132)  with  the  m-1  polynomials  in  its  denominator  is  required. 
The  numerical  poles/zeros  cancellations  are  not  exact  due  to  the  inevitable  computer  round-off.  The  m-1 
such  cancellations  of  each  zero  of  the  numerator  with  the  m-1  poles  of  the  denominator  must  occur.  Thus, 
in  order  to  recover  numerical  accuracy  one  must  factor  all  the  numerator  and  denominator  polynomials  of 
P  and  then  check  out  the  inevitable  inexact  cancellations  as  demonstrated  in  Fig.  56. 
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Roots  of  High  Precision  Polynomial  (Mathematica). 
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Roots  of  Standard  Precision  Polynomial  (MATRIXx,  Control-C,  MATLAB,  Macsyma). 

Figure  56:  Figures  by  Sating,  (1992). 


5.0  APPLICATIONS 
5.1  MIMO  QFT  Example 

Design  Problem  —  A  2x2  analog  flight  control  system  (see  boxed  in  loops  of  Fig.  40)  illustrates  the 
increased  accuracy  and  the  efficiency  achieved  by  use  of  the  MIMO  QFT  CAD  package  developed  by  Sating, 
(1992).  This  CAD  package  is  a  straight  forward  method  for  designing  an  analog  or  discrete  MIMO  control 
system.  The  original  design  was  done  by  Arnold,  (1984)  but  was  redesigned  by  Sating  utilizing  his  CAD 
package  [Sating,  (1992)].  The  specifications  are:  to  design  a  robust  analog  aircraft  flight  control  system, 
provide  the  required  stability,  and  to  satisfy  the  time  domain  performance  requirements.  Four  flight 
conditions  (Table  4)  and  six  aircraft  failure  modes  (Table  5)  are  specified.  Table  6  lists  the  resulting  set  of 
24  plant  cases  that  incorporate  these  flight  conditions  and  failure  modes. 
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Table  4:  Flight  Conditions 


Flight  Condition 

Aircraft  Parameters 

Mach 

Altitude 

1 

0.2 

30 

2 

0.6 

30,000 

3 

0.9 

20,000 

4 

1.6 

30,000 

Failure  Mode 

Failure  Condition 

1 

Healthy  aircraft 

2 

One  horizontal  tail  fails 

3 

One  flaperon  fails 

4 

One  horizontal  tail  and  one 
flaperon  fail:  same  side 

5 

One  horizontal  tail  and  one 
flaperon  fail:  opposite  sides 

6 

Both  flaperons  fail 

Table  6:  Plant  Models 


Failure 

Flight  Condition 

Mode 

1 

2 

3 

4 

1 

#1 

#7 

#13 

#19 

2 

#2 

#8 

#14 

#20 

3 

#3 

#9 

#15 

#21 

4 

#4 

#10 

#16 

#22 

5 

#5 

#11 

#17 

#23 

6 

#6 

#12 

#18 

#24 

The  P.S.  are  to  design  for  a  stability  of  y  =  45°  phase  margin  for  each  of  the  two  feedback  loops. 
The  frequency  domain  perfonnance  specifications,  when  met,  result  in  the  desired  closed  loop  system 
performance  in  the  time  domain  and  are  shown  as  dashed  lines  on  the  Bode  plots  of  Fig.  57. 


Figure  57:  The  2x2  System  Performance: 
the  Specified  Performance  Bounds  in  Dashes  and  tjj  Responses. 
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Design  Process  —  The  specifications,  plant  models  for  the  24  cases,  and  the  weighting  matrix,  W, 
are  entered  into  the  QFT  CAD.  If  the  plant  matrix  P  is  not  square  then  an  effective  square  plant  matrix 
Pe  =  PW  is  obtained  in  order  to  obtain  a  QFT  design.  Automated  features  accessed  through  the  designer 
CAD  interface  result  in  synthesizing  the  compensators  g  ]  (s)  and  g2(s).  The  nominal  loop  transmission 
functions  Ll0(s)  =  g  (s)qllo(s)  and  L2o(s)  =  g2(s)q22o(s)  are  synthesized  (or  shaped)  so  that  they  satisfy  their 
respective  stability  bounds  and  their  respective  optimal  bounds  Blo(jeoi)  and  B2o(jcoj).  Note  that  qUo  and 
q2,  are  the  nominal  plant  transfer  functions. 

Validation  Check  —  A  first  step  in  the  validation  process  is  to  obtain  plots  of  the  loop  transmission 
functions  L2l(s),  where  i  =  1,  ...  ,  24,  for  all  24  cases  on  the  N.C.  These  plots  are  obtained  by  a  CAD  routine, 
see  Fig.  58,  for  the  purpose  of  a  stability  check.  Note,  that  all  the  plants  are  minimum-phase.  None  of  the 
cases  violate  the  ML  stability  contour  to  accept  the  consequences  of  violating  the  disturbance  bound  for 
co  =  2  r/s.  With  Llo(s)  and  L2o(s)  synthesized,  the  CAD  package's  automated  features  expedite  the  design  of 
the  prefilters  fn(s)  and  f22(s). 


Figure  58:  NC  Stability  Validation  Check  of  the  l_2  (j<o)  Plots:  24  Open-loop  Transmissions. 
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As  a  second  step  in  the  design  validation  process  the  2x2  array  of  Bode  plots,  shown  in  Fig.  57,  are  generated 
showing  on  each  plot  the  24  possible  closed  loop  transmissions  from  an  input  to  an  output  of  the  completed 
designed  system.  The  consequence  of  violating  the  channel  2  disturbance  bound,  for  co  =  2  r/s,  is  seen  where 
the  closed  loop  transmissions  violate  the  bound  b2i.  This  bound  is  denoted  by  a  dashed  line.  Beginning  at 
about  co  =  2  r/s  a  violation  of  the  performance  bounds  during  loop  shaping  may  result  in  violation  of  the 
performance  specifications  for  the  closed  loop  system.  There  is  a  violation  of  the  bn  bound  of  Channel  1. 
The  penetration  of  the  upper  tracking  bound  at  co;  «  0.1  r/s  is  due  to  the  phugoid  mode  which  is  not  a 
problem;  i.e.,  does  not  affect  the  P.S. 

Final  Check  —  Figure  57  shows  a  robust  design  has  been  achieved  for  this  2x2  MIMO  analog  flight  control 
system.  The  time  domain  results,  although  not  drawn,  meet  all  performance  specifications. 

5.2  Other  Examples 

The  following  is  a  list  of  other  QFT  designed  examples  that  are  mentioned  by  Floupis,  et  al.  (1999): 
QFT  And  Robust  Process  Control,  Idle  Speed  Control  For  Automotive  Fuel  Injected  Engine, 
Welding  Control  Systems,  Control  System  For  An  Actuator  Plant,  VISTA  F-16  Flight  Control  System 
(Including  Configuration  Variation),  Design  Of  Flight  Control  Laws  For  Aircraft  With  Flexible  Wings 
Using  Quantitative  Feedback  Theory,  Robot  Controllers,  and  Wastewater  Treatment  Control  System. 


6.0  SUMMARY 

The  engineers  are  applying  the  results  of  the  scientific  method  to  achieving  solutions  for  the  2 1  st  century 
real  world  problems.  The  QFT  approach  has  the  advantage  that  it  is  close  to  the  engineers’  existing 
experience  in  conventional  (classical)  design  methods.  It  also  provides  the  facilities  to  deal  with  structured 
and  unstructured  parametric  uncertainty  which  appear  not  to  be  available  in  traditional  methods.  QFT  has 
the  attractive  features  of  providing  a  link  with  existing  techniques  whilst  at  the  same  time  providing  many 
of  the  advanced  features  that  are  needed  for  the  2 1st  century  high  performance  systems. 

The  QFT  Control  System  Design  Process  typifies  the  “Bridging  the  Gap”  (Houpis,  2002).  As  an  example, 
a  real-world  QFT  robust  multivariable  control  system  design  problem  for  an  Unmanned  Research  Vehicle 
(URV)  is  presented  by  Houpis,  et  al.  (1999)  which  involves  the  elements  in  “Bridging  the  Gap” 
as  indicated  in  Fig.  59.  Note  that  at  each  step  of  the  design  process,  as  denoted  by  the  rectangular  boxes, 
the  engineer  can  ascertain  if  the  P.S.  are  still  being  satisfied.  If  not,  as  indicated  by  the  curve  line 
emanating  from  each  box  that  terminates  on  the  “Redesign”  box  the  engineer  does  not  have  to  complete 
the  entire  design  process  before  needing  to  do  a  “redesign.” 
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ORGANIZATION 


Engineering  Visualization 

Used  to  better  understand  results 
of  simulations  and  System  tests 


r 

Engineering  Interactive  Simulation 

•User  supplies  commands  and  then  can 
react  to  resulting  dynamic  behavior 
•Gives  a  better  understanding  of  control 
system  operation 


Hardware-in-the-Loop 

Simulation/Implementation 

•Real-time  operation  of  control  algorithm 
•Noise  corrupted  measurements  available 
for  feedback 

•Computation  cycle  time/Sampling  Rate 
•Quantization  Error,  Warping 


Figure  59:  The  QFT  Control  System  Design  Process:  Bridging  the  Gap. 
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